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Preface 

The  preprint  entitled  "Laser  Optogalvanic  and  Fluorescence  Studies  of  the 
Cathode  Region"  forms  the  body  of  this  report.  The  paper  summarizes  more  than 
three  years  of  research  on  the  cathode  fall  region  of  glow  or  diffuse 
discharges.  Although  complete  self-consistent  models  of  the  cathode  region  are 
not  yet  available,  we  have  made  significant  progress  toward  this  goal.  The 
experimental  and  theoretical  results  in  the  preprint  provide  a  detailed 
miscrosopic  picture  of  the  cathode  region. 
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Abstract 

Various  laser  diagnostics  are  used  to  study  the  cathode  fall  and  negative 
glow  regions  of  a  He  glow  discharge  with  a  cold  A1  cathode.  The  electric  field 
and  absolute  metastable  densities  are  mapped  and  the  gas  temperature  is 
measured  over  a  range  of  current  densities  from  a  near  normal  (173 V)  to  a 
highly  abnormal  (600V)  cathode  fall.  These  measurements  are  analyzed  to  yield 
the  current  balance  at  the  cathode  surface,  the  ionization  rate  in  the  cathode 
fall  region,  and  the  metastable  production  rate  in  the  cathode  fall  and 
negative  glow  regions.  The  experimental  results  compare  favorably  with  the 
results  of  Monte  Carlo  simulations.  The  density  and  temperature  of  the  low 
energy  electron  gas  in  the  negative  glow  is  determined  by  combining  information 


from  the  experiments  and  Monte  Carlo  simulations. 


-2- 


I.  Introduction 

There  has  in  recent  years  been  a  very  substantial  effort  toward  the 
development  of  a  quantitative  microscopic  understanding  of  the  cathode  region 
of  glow  discharges.1  ®  This  region  is  of  great  practical  interest  because  of 
its  importance  in  plasma  processing  and  in  pulsed  power  devices.  Plasma 
processing  includes  applications  of  glow  discharges  in  ion  etching,  thin  film 
deposition,  and  plasma  treating  of  surfaces.  The  stability  of  pulsed  gas 
discharge  lasers  and  other  pulsed  power  devices  is  largely  determined  by  the 
stability  of  the  cathode  fall  region. 

The  cathode  region  is  also  of  great  fundamental  interest.  The 
hydrodynamic  approximation,  which  is  often  very  useful  in  the  central  region  of 
a  discharge,  fails  in  the  cathode  region.  If  the  local  E/N  (electric  field  to 
gas  density  ratio)  determines  the  velocity  dependence  of  the  electron 
distribution  function,  then  the  electrons  are  in  hydrodynamic  equilibrium.^ 
The  lack  of  hydrodynamic  equilibrium  in  the  cathode  fall  region  is  caused  by 
the  large  and  rapidly  changing  E/N,  and  by  the  proximity  of  the  electrode.  The 
cathode  region  includes  both  the  cathode  fall  and  negative  glow  regions.  The 
negative  glow  has  a  high  density  of  low  energy  electrons  and  very  small 
electric  fields.  It  also  is  a  nonhydrodynamic  region  because  of  a  low  density 
of  high  energy  "beam"  electrons  injected  from  the  cathode  fall  region. 
Unfortunately  the  extensive  data  base  of  electron  transport  coefficients  such 
as  drift  velocities,  diffusion  coefficients,  and  Townsend  coefficients  is  not 
directly  applicable  in  modelling  the  cathode  region.^  These  coefficients  are 


generally  measured  in  drift  tube  experiments  in  which  considerable  effort  is 
devoted  to  insuring  that  the  hydrodynamic  approximation  holds. 
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Nonequilibrium  regions  are  receiving  increased  attention  from  theorist  and 
experimentalist.  Because  of  their  nonequilibrium  nature,  electrons  in  the 
cathode  fall  region  are  often  modelled  using  either  a  full  Boltzmann  equation 
analysis  or  a  Monte  Carlo  simulation,  A  somewhat  less  realistic  but  much 
faster  approach  is  based  on  a  model  distribution  function  with  position 

1  O 

dependent  parameters.  The  spatial  dependences  of  the  parameters  such  as 
electron  density,  average  velocity,  and  average  energy,  are  determined  by 
solving  moments  of  the  Boltzmann  equation  which  are  coupled  ordinary 
differential  equations.  These  approaches  might  in  general  be  called 
nonequilibrium  fluid  models.  Single  beam  and  multiple  beam  models  are  examples 

1  ll  1C 

of  this  general  approach.  ’  J  This  approach  is  attractive  because  it  is 
usually  easier  and  faster  to  solve  moments  of  the  Boltzmann  equation  which  are 
coupled  ordinary  differential  equations  than  it  is  to  solve  the  full  Boltzmann 
equation  which  is  a  partial  differentio-integral  equation.  The  nonequilibrium 
fluid  approach  can  in  principal  be  made  quite  realistic  if  one  is  sufficiently 
clever  in  the  design  of  a  model  distribution  function  for  the  electrons. 

A  more  complete  solution  will  require  a  calculation  of  self  consistent 
space  charge  electric  fields.  This  will  involve  coupling  Poisson's  equation 
and  one  (or  more)  moment  equation(s)  for  ions  to  the  moment  equations  for 
electrons.  Consider  a  model  which  includes  three  moments  of  the  Boltzmann 
equation  for  electrons,  Poisson's  equation,  and  a  single  continuity  equation 
for  ions.  This  model  requires  the  simultaneous  solution  of  five  coupled  first 
order  differential  equations.  Four  natural  boundary  conditions  include 
specifying  the  electron  distribution  function  at  the  cathode  and  specifying  an 
electron  emission  coefficient  for  ion  bombardment  of  the  cathode.  The 
"missing"  condition  is  unfortunately  not  a  boundary  condition  but  rather  an 


extremum  condition.  The  physically  correct  solution  will  maximize  the  current 
at  fixed  voltage  or  minimize  the  voltage  at  fixed  current. 

The  problem  described  in  the  preceeding  paragraph  is  obviously  quite 
difficult.  It  is  also  incomplete  because  of  the  important  role  played  by 
metastable  atoms,  and  UV  or  VUV  photons  in  releasing  electrons  from  the 
cathode.  Appropriate  balance  equations  for  these  neutral  particles  must  also 
be  included  in  the  problem. 

A  complete  solution  to  the  cathode  fall  problem  is  beyond  the  scope  of 
this  paper.  This  paper  describes  extensive  experimental  results  from  recently 
developed  laser  optogalvanic  diagnostics,  and  more  traditional  lasei — induced 
fluorescence  and  laser  absorption  diagnostics.1^  These  experimental  results 
include  electric  field  maps,  gas  temperature  (density)  measurements,  and 
absolute  metastable  density  maps  in  a  He  glow  discharge  over  a  range  of  current 
densities  from  a  near  normal  cathode  fall  of  173V  to  a  highly  abnormal  cathode 
fall  of  600V.  These  experimental  results  are  analyzed  in  detail  to  provide 
information  on  electron  impact  ionization  and  excitation  in  the  cathode  fall 
and  negative  glow  regions.  Monte  Carlo  simulations  of  electron  avalanches  in 
the  cathode  region  are  compared  to  the  experimental  results.  The  Monte  Carlo 
code  is  based  on  Boeuf  and  Marode's  variation  of  the  null  collision  technique 
for  non-uniform  fields.1^  The  agreement  of  the  Monte  Carlo  results  and  the 
experimentally  derived  ionization  and  excitation  rates  is  generally  good. 
Finally  the  experimental  results  are  combined  with  Monte  Carlo  simulations  to 
infer  the  density  and  temperature  of  the  low  energy  electron  gas  in  the 
negative  glow.  This  unexpected  information  is  largely  based  on  observation  of 
a  suppression  of  the  2^  metastable  density  in  the  negative  glow  due  to 
metastable  spin  conversion  by  low  energy  electron  collisions. 
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The  results  of  this  paper  can  be  viewed  in  several  ways:  (1)  as  a  test  of 
our  understanding  of  the  key  physical  processes  in  the  cathode  region,  (2)  as  a 
set  of  benchmark  experiments  for  further  modelling  of  the  cathode  region,  or 
(3)  as  a  starting  point  for  more  ambitious  experiments  to  study  electron 
emission  from  the  cathode.  We  plan  to  make  in-situ  determinations  of  the 
relative  importance  of,  and  coefficients  for,  electron  emission  due  to  ion, 
metastable  atom,  and  UV  or  VUV  photon  bombardment  of  the  cathode. 

The  structure  of  this  paper  is  as  follows.  Section  II  is  a  description  of 
the  experimental  apparatus.  Sections  III  and  IV  are,  respectively,  a 
presentation  of  the  electric  field  and  gas  density  measurements,  and  an 
analysis  of  this  data  to  yield  information  on  ionization  in  the  cathode  fall 
region.  The  field  maps  and  gas  density  measurements  are  combined  with  an 
analytic  treatment  of  ion  transport  in  the  cathode  fall  region  to  determine  the 
current  balance,  or  ratio  of  ion  to  electron  current  at  the  surface  of  the 
cathode.  This  ratio  is  found  to  be  3.3  and  is  independent  of  the  total 
discharge  current.  Section  V  and  VI  are  respectively,  a  presentation  of  the 
metastable  density  maps  and  an  analysis  of  metastable  diffusion  and  kinetics  in 
the  cathode  fall  and  negative  glow  regions.  The  analysis  of  the  data  yields 
?1S  and  2^S  metastable  production  rates  and  a  rate  for  2^  metastable  quenching 
due  to  low  energy  electron  collisions  in  the  negative  glow.  Sections  VII  and 
VIII  are,  respectively,  a  discussion  of  the  Monte  Carlo  code  including  electron 
impact  cross  sections  and  a  comparison  of  the  Monte  Carlo  results  with  the 
experimental  results.  Section  VIII  also  includes  an  analysis  of  the  density 
and  temperature  of  the  low  energy  electrons  in  the  negative  glow.  Finally 
Section  IX,  includes  a  summary,  some  conclusions,  as  well  as  speculations  on 


future  work. 


II.  Experimental  Apparatus 


Figure  i  is  a  schematic  of  the  experimental  apparatus.  The  discharge  used 
in  these  studies  is  produced  between  flat  circular  A1  electrodes  3*2  cm  in 
diameter  and  separated  by  0.62  cm.  The  electrodes  are  water-cooled  to  minimize 
gas  heating.  The  discharge  tube  is  made  primarily  of  glass  and  stainless 
steel.  Most  of  the  large  seals  are  made  with  knife  edge  flanges  on  Cu  gaskets. 
The  only  exceptions  are  the  high  vacuum  epoxy  seals  around  the  fused  silica 
brewster  windows.  A  liquid  N2  trapped  diffusion  pump  evacuates  the  tube  to 

_Q 

2x10  Torr.  When  no  liquid  N2  is  in  the  trap  an  ion  pump  maintains  the  vacuum 
to  prevent  oil  from  back-diffusing  into  the  system.  The  leak  rate  into  the 
sealed  discharge  tube  is  approximately  3x10_1<  Torr /day.  For  discharge 
operation  ultra  high  purity  (0.999999)  He  is  slowly  flowed  through  the  system. 

A  capacitive  manometer  monitors  the  pressure  which  is  maintained  at  3.5  Torr. 
The  He  first  passes  through  a  cataphoresis  discharge  to  remove  any  residual 
contaminants  before  entering  the  main  discharge  tube.  Emission  spectra  reveal 
only  weak  A1  and  H  impurity  lines.  Although  He  is  much  less  effective  at 
sputtering  than  heavy  inert  gases,  some  sputtering  of  the  A1  cathode  does 
occur.  We  suspect  that  the  slow  erosion  of  the  cathode  gives  rise  to  both  the 
weak  A1  and  H  impurity  lines.  Many  metals  absorb  hydrogen.  The  surface  of  the 
A 1  cathode  is  cleaned  in-situ  by  running  an  Ar  discharge.  Argon  is  very 
effective  at  sputtering  most  metals  including  Al.  A  carefully  prepared  cathode 
provides  discharge  V-I  characteristics  which  are  stable  to  a  few  percent  over  a 
period  of  a  month. 

Experiments  are  carried  out  over  a  range  of  discharge  current  densities 
from  O.190  mA/cm  to  1.50  mA/cm  .  The  low  current  density  corresponds  to  a 
near  normal  cathode  fall  voltage  of  173  V  and  the  high  current  density 
corresponds  to  a  highly  abnormal  cathode  fall  voltage  of  600  V.  The  discharge 


current  is  spread  uniformly  across  the  surface  of  the  electrodes.  The  absence 
of  significant  fringing  was  verified  by  segmenting  the  3*2  cm  diameter  cathode 
into  a  1 .6  cm  diameter  disc  and  a  close  fitting  annulus  with  an  outside 
diameter  of  3.2  cm.  Each  part  of  the  cathode  was  maintained  at  the  same 
potential  during  operation  and  the  average  current  density  on  each  was 
measured.  The  current  density  was  found  to  be  uniform  across  the  cathode. 

This  measurement  is  one  of  the  justifications  for  using  one  dimensional  models 
of  the  cathode  fall  region.  The  segmented  cathode  was  replaced  by  a  solid 
cathode  after  the  current  density  measurements. 

The  laser  used  is  a  N2  laser  pumped  dye  laser.  Several  different 
configurations  for  the  dye  laser  are  used,  depending  on  the  specific 
experiment.  The  dye  laser  is  tunable  over  a  wavelength  range  360-700  nm. 
Frequency  doubling  crystals  are  used  to  extend  this  range  down  to  200  nm.  For 
electric  field  measurements  the  laser  is  frequency  doubled  in  order  to  excite 
Rydberg  levels  from  the  heavily  populated  metastable  levels.  The  dye  laser  is 
operated  with  a  bandwidth  of  0.2  to  0.H  cm-1.  The  spectral  widths  of 
components  of  Rydberg  atom  Stark  manifolds  are  often  much  broader  than  a  few 
tenths  of  a  cm  1  due  to  the  spatial  gradients  of  the  electric  field  in  the 
cathode  fall.  We  find  it  unnecessary  to  reduce  the  dye  laser  bandwidth  below 
0.2  cm  1  for  electric  field  measurements  using  Rydberg  atoms.  The  dye  laser 
bandwidth  is  reduced  to  0.01  cm-1  (300  MHz)  for  measuring  the  Floppier  width  of 
a  transition  between  non-Rydberg  levels.  The  dye  laser  is  operated  with  a 
bandwidth  of  0.2  cm  ^  to  0.M  cm  ^  for  laser  induced  fluorescence  studies  of  the 
spatial  distribution  of  metastable  atoms.  The  absolute  metastable  density  is 
provided  by  an  absorption  experiment.  The  pulsed  dye  laser  power  must  be 
substantially  reduced  to  avoid  saturation  in  the  absorption  experiments.  We 
also  found  it  convenient  to  reduce  the  dye  laser  bandwidth  to  0.017  cm-1.  A 
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very  small  instrumental  contribution  simplifies  analysis  of  the  absorption 
data . 

Figure  1  shows  all  three  detection  schemes  used  in  these  experiments. 
Optogalvanic  detection  involves  using  the  boxcar  averager  to  monitor  the 
current  through  the  discharge  by  observing  voltage  changes  across  the  263  ft 
resistor  between  the  cathode  and  ground.  This  detection  scheme  is 
particularity  sensitive  in  the  cathode  fall  region  because  optogalvanic  effects 
are  naturally  amplified  by  the  cathode  fall.18  Optogalvanic  detection  is 
essential  for  observing  transitions  to  Rydberg  levels.  The  fragile  Rydberg 
atoms  do  not  live  long  enough  in  order  to  fluoresce  in  the  visible  or  UV;  they 
are  rapidly  coll isionally  ionized  through  associative  ionization  and  other 
mechanisms.  The  oscillator  strengths  for  transitions  from  the  heavily 
populated  metastable  levels  to  the  Rydberg  levels  are  very  small  and  a  direct 
observation  of  absorption  would  be  extremely  difficult.  Optogalvanic  detection 
is  also  used  in  the  Doppler  width  measurements  reported  in  this  work. 
Optogalvanic  detection  is  not  suited  to  map  the  metastable  density.  The 
efficiency  of  the  optogalvanic  effect  is  strongly  position  dependent  in  the 
cathode  region.  Laser-induced  fluorescence  provides  the  most  convenient 
detection  method  for  mapping  the  spatial  distribution  of  meta3tables.  The 
fluorescence  is  detected  with  a  photomultiplier  after  passing  through  an 
interference  filter  which  blocks  much  of  the  spontaneous  emission  from  the 
discharge  and  much  of  the  laser  light  scattered  from  windows,  etc.  The  least 
sensitive  detection  method  used  in  the  experiments  is  direct  absorption.  This 
method  does  provide  a  convenient  absolute  calibration  for  the  metastable 
density  maps.  The  photodiode  shown  in  Fig.  1  is  used  in  absorption 


experiments . 


-9- 


In  order  to  perform  spatially  resolved  measurements  without  disturbing 
laser  alignment  the  discharge  is  mounted  on  a  precision  translation  state. 

Good  spatial  resolution  is  achieved,  when  necessary,  by  passing  the  laser  beam 
through  a  narrow  slit  and  imaging  the  slit  into  the  discharge  region. 

III.  Electric  Field  and  Gas  Density  Measurements 
Space  charge  electric  fields  in  gas  discharges  have  traditionally  been 
measured  using  electron  beam  deflection.  A  laser  induced  fluorescence 
technique  for  field  measurements  was  recently  developed  by  Moore,  Davis,  and 
Gottseho.1^  a  laser  technique  based  on  optogalvanic  detection  of  Rydberg  atoms 
was  recently  developed  by  Doughty,  Salih  and  Lawler.1^*'20  Both  laser  techniques 
have  a  number  of  important  advantages  over  the  traditional  electron  beam 
deflection  technique.  The  laser  techniques  are  useful  at  higher  pressures  and 
discharge  current  densities.  It  is  easier  to  use  a  spectroscopic  technique  in 
a  high  purity  discharge  system.  The  laser  techniques  are  truly  nonperturbing 
when  performed  with  a  nanosecond  pulsed  laser  because  the  field  is  measured 
when  the  atoms  or  molecules  absorb  the  laser  light,  and  any  perturbation  to  the 
discharge  fields  occurs  on  a  longer  time  scale.  Lasers  have  the  potential  for 
making  measurements  with  a  few  microns  spatial  resolution  and  nanosecond 
temporal  resolution.  The  technique  based  on  Rydberg  atoms  has  broad 
applicability  because  all  atoms  and  molecules  have  Rydberg  levels.  This 
technique  also  has  a  wide  dynamic  range  because  one  can  choose  an  appropriate 
principal  quantum  number.  The  Rydberg  atom  technique  has  been  used  to  measure 
fields  as  small  as  10  V/cm  ±  1  V/cm.22 

The  technique  based  on  optogalvanic  detection  of  Rydberg  atoms  excited 
with  a  single  laser  is  used  in  this  investigation.1^  A  technique  involving 
Intersecting  laser  beams  for  pinpoint  field  measurements  has  also  been 
demonstrated.  The  discharge  studied  in  this  experiment  is  sufficiently  one 


dimensional  that  the  simpler  single  laser  technique  is  adequate.  In  addition 
to  the  advantages  described  in  the  preceeding  paragraph,  a  Rydberg  atom 
technique  offers  an  advantage  in  ease  of  interpretation.  The  experimental 
Stark  spectra  are  analyzed  by  comparing  to  theoretical  Stark  maps. 
Straightforward,  ab  initio  calculations  of  these  Stark  maps  to  accuracy  are 
possible  for  many  atoms.  An  excellent  procedure  for  calculating  the  Stark  maps 
is  described  in  detail  by  Zimmerman,  Littman,  Kash,  and  Kleppner.2^  The 
procedure  used  in  this  work  is  a  simplified  version  of  that  described  by 
Zimmerman  et .  al.  The  main  simplification  is  that  we  use  corrected  hydrogenic 
radial  matrix  elements  instead  of  performing  a  numerical  integration  for  radial 
matrix  elements. 

A  Rydberg  atom  is  an  atom  with  one  of  its  electrons  excited  to  a  high 
lying  energy  level  characterized  by  a  large  principal  quantum  number,  n.  As  n 
increases  many  atomic  properties  tend  towards  those  of  hydrogen.  The  energy  of 
Rydberg  states  is  given  by  the  empirical  formula. 


V  ■  “ion- 


(n-5^)' 


(1) 


where  Wion  is  the  ionization  potential,  R0  is  the  Rydberg  constant,  l  is  the 
angular  momentum  quantum  number,  and  6^  is  the  quantum  defect.  This  expression 
describes  hydrogenic  energy  levels  for  6^  -  0.  The  quantum  defect  reflects  the 
degree  of  interaction  of  the  outer  valence  electron  with  the  core  electrons. 

For  large  values  of  l  the  valance  electron’s  orbit  is  nearly  circular,  so  it 
does  not  interact  to  any  significant  degree  with  the  core.  Consequently,  this 
electron  sees  essentially  a  Coulomb  potential  implying  6^  is  vanishingly  small. 
Low  angular  momentum  orbits,  such  a3  s  and  p  orbits,  do  penetrate  the  core, 
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thus  interacting  strongly  with  the  charges  there.  On  average,  an  s  orbit  will 
see  a  net  core  charge  greater  than  one  because  of  its  penetration,  so  this 
level  will  be  more  tightly  bound  than  the  other  levels  of  the  same  principal 
quantum  number  and  the  quantum  defect  will  be  larger.  These  qualitative 
descriptions  explain  why  6^  is  weakly  dependent  on  n  but  strongly  dependent  on 
£. 

Mixing  of  states  with  different  principal  quantum  numbers  occurs  if  both 
the  field  and  principal  quantum  number  are  sufficiently  large.  This  is  an 
undesirable  complication  if  one  is  studying  a  discharge  and  simply  using 
Rydberg  atoms  as  a  probe.  It  can  usually  be  avoided  by  choosing  an 
appropriate,  somewhat  lower  principal  number  for  very  high  fields.  In  heavy 
atoms  with  large  quantum  defects  it  may  be  necessary  to  include  states,  which 
are  nearby  at  zero  field,  but  have  different  principal  quantum  numbers.  In  He 
all  quantum  defects  are  less  than  0.5. 

The  energy  perturbations  are  calculated  by  diagonalizing  an  (n-m)  by  (n~m) 

Hamiltonian  matrix,  where  m  is  the  component  of  angular  momentum  along  the 

field  axis.  The  basis  set  is  made  of  eigenfunctions  of  the  unperturbed 

Hamiltonian.  The  diagonal  matrix  elements  are  the  unperturbed  energies  of  the 

levels  as  given  in  Eq.  1.  The  quantum  defect,  6%,  can  be  deduced  from  Moore's 

pH 

tables  of  energy  levels  .  The  off-diagonal  elements  vanish  except  those 
connecting  states  for  which  i  differs  by  one.  These  matrix  elements  are 
approximately  the  hydrogenic  matrix  elements  given  by  the  expression: 

<n  &  m|eE«r|n  £-1  m> 

-  (3/Z)eEa0n[ (nZ-»2Ht2-»l)]1/2. 

H£2-1 


(2) 


whore  e  is  the  unit  charge,  and  aQ  is  the  Bohr  radius.  The  dipole  matrix 
elements  above  must  be  corrected  to  account  for  core  penetration  effects.  A 
correction  factor  of  order  one  has  been  tabulated  as  a  function  of  quantum 
defect  by  Edmonds,  et.  al.^ 

Figure  2  presents  the  results  of  applying  the  procedure  outlined  above  to 

the  n  =  11,  m  =  0,  singlet  levels  of  He.  The  diagonal  elements  are  determined 

from  the  splittings  of  the  I^S,  1 1 1 P ,  I^D,  1 1 1 F ,  and  l^G  levels,  whereas  the 

H,  I,  K,  L,  M  and  N  levels  are  approximated  as  being  degenerate  with  the  l^G 
2H 

level.  A  5%  correction  is  applied  to  the  S-P  matrix  elements  to  account  for 

or 

core  penetration.  J  Corrections  to  the  other  elements  are  generally  less  than 
1£.  Figure  2a  is  a  plot  of  the  energies  of  the  m  =  0  states  for  the  n  =  11 
levels  as  a  function  of  applied  electric  field.  The  energies  are  the 
eigenvalues  of  the  diagonalized  Hamiltonian.  Only  ten  components  appear  in 
this  figure  since  the  I^S  level  has  a  large  quantum  defect  and  is  23  cm  1 
below  this  pattern.  This  level  does  not  mix  strongly  with  the  other  n  =  11 
levels  at  the  fields  indicated  in  the  figure.  For  fields  above  200  V/cm  the 
pattern  clearly  reflects  a  linear  Stark  effect. 

Figure  2b  is  a  plot  of  the  expected  intensity  of  selected  Stark  components 
from  Fig.  2a  as  a  function  of  electric  field.  Suppose  transitions  between  a 
low-lying  non-Rydberg  ^  level  and  the  n  »  11  Rydberg  states  are  driven.  With 
no  externa]  electric  field,  only  the  transition  will  be  observed.  The 

effect  of  the  field  is  to  mix  a  fraction  of  the  1 1 1 P  state  into  each  n  -  11 
state.  This  fraction  Is  determined  by  squaring  the  appropriate  element  from  an 
eigenvector  of  the  diagonalized  Hamiltonian.  The  fraction  is  a  relative 
linestrength  for  a  transition  to  that  component  of  the  Stark  manifold.  The 
relative  linestrengths  of  the  components  labeled  1 ,  3.  and  8  are  plotted  in 


Fig.  ^b.  Such  information  is  useful  in  determining  electric  fields  in  the  ease 
where  the  splittings  between  individual  components  are  small. 

In  order  to  use  Rydberg  atoms  to  measure  discharge  electric  fields  these 
states  must  first  be  produced  (there  is  no  appreciable  steady  state  population 
in  the  cathode  fall  region)  and  then  detected.  The  narrow  bandwidth  and  high 
power  characteristics  of  the  laser  make  it  an  ideal  tool  to  excite  atoms  from  a 
well  populated  lower  level  to  a  specific  Rydberg  level.  Since  the  oscillator 
strength  decreases  as  1/n^  it  is  very  difficult  to  directly  detect  the  weak 
absorption.  Laser  induced  fluorescence  from  these  atoms  is  also  difficult  to 
detect,  because  in  the  1  -  10  Torr  pressure  regime  the  Rydberg  atoms  do  not 
survive  long  enough  to  radiate.  They  are  collisionally  ionized  via  collisions 
with  ground  state  atoms.  Associative  ionization: 

He* ( nil > 3s)  +  HeO^)  -  He!,  ♦  e"  (3) 

is  likely  important  for  Rydberg  levels  below  the  ionization  limit  by  more  than 
thermal  energies.  The  absorption  of  laser  radiation  produces  excess  ionization 
in  the  discharge  and  thus  produces  an  optogalvanic  effect.  An  optogalvanic 
effect  is  detectable  as  a  perturbation  of  the  discharge  current  and  can  in 
general  be  either  positive  (increased  current)  or  negative.  Absorption  of 
radiation  may  in  some  cases  cause  a  reduction  of  ionization  in  the  discharge 
and  a  decrease  in  discharge  current.  Optogalvanic  spectroscopy  is  uniquely 
suited  to  the  detection  of  Rydberg  atoms  in  a  discharge  environment. 

The  cathode  fall  region  is  particularly  well  suited  to  optogalvanic 
detection.  If  an  extra  ion-electron  pair  is  introduced  into  the  cathode  fall 
region  the  electron  is  accelerated  by  the  field,  thus  producing  additional 
ionization.  This  amplification  can  result  in  optogalvanic  effects  in  the 


cathode  fall  region  that  are  100  times  more  sensitive  than  optogalvanic  effects 

1  ft 

in  the  positive  column. 

If  the  lower  level  from  which  the  transition  is  driven  is  metastable,  the 
ionization  process  results  in  the  destruction  of  the  metastable  atom.  The 
destruction  of  atoms  in  these  levels  can  produce  a  negative  optogalvanic  effect 
since  metastables  diffusing  to  the  cathode  have  a  large  probability  of 
releasing  electrons  from  the  cathode.  It  would  seem  that  this  negative 
optogalvanic  effect  might  cancel  the  positive  optogalvanic  effect  which  is  due 
to  the  extra  electron  produced  by  associative  ionization.  The  time  scales  for 
the  two  signals,  however,  are  vastly  different.  The  diffusion  time  of  the 
metastables  to  the  cathode  can  be  several  hundred  microseconds,  whereas 
associative  ionization  produces  a  signal  in  less  than  1  microsecond.  When 
produced  by  a  pulsed  laser  the  positive  portion  of  the  transient  optogalvanic 
signal  is  well  isolated  and  can  be  easily  monitored  with  a  boxcar  averager. 

The  critical  consideration  in  measuring  discharge  electric  fields  is 
choosing  an  appropriate  transition  to  produce  the  Rydberg  atoms.  For  noble 
gases  in  general  and  He  in  particular,  transitions  from  the  ground  state  to 
upper  levels  are  inaccessible  with  present  dye  lasers.  These  gases,  however, 

have  metastable  levels  which  are  well  populated  in  the  discharge.  In  He,  both 

10  11 
<  S  and  P-’S  levels  are  metastable.  The  2  S  to  11  P  transition  at  321  nm  is  the 

most  suitable  under  the  conditions  of  this  experiment.  Although  the  ?  level 

has  a  higher  population  than  the  21S,  the  n^P  levels  have  a  much  higher  quantum 

defect  than  the  n'p's  (0.068  and  -0.012  respectively)24.  The  larger  quantum 

defect  necessitates  exciting  to  a  level  with  a  larger  n  to  observe  a  linear 

Stark  effect  for  a  given  field.  The  oscillator  strength  goes  as  1/n^  while  the 

number  of  components  increases  as  n,  so  the  amount  of  signal  in  each  component 


Radiation  from  the  dye  laser  is  frequency  doubled  using  a  potassium 
d i hydrogen  phosphate  crystal.  The  frequency  doubled  beam  is  polarized  normal 
to  the  surface  of  the  cathode  so  that  only  Am  =  0  transitions  are  excited. 

This  results  in  only  m  =  0  Rydberg  states  because  the  lower  level  is  a  1S.  The 
beam  is  focussed  in  the  center  of  the  discharge  to  a  strip  approximately 
0.01  cm  wide  and  1  cm  long  parallel  to  the  surface  of  the  electrodes.  As  the 
laser  is  scanned  through  the  Stark  manifold,  the  current  is  monitored  across  a 
263 s;  resistor  between  the  cathode  and  ground.  The  optogalvanic  signal  is 
processed  using  a  boxcar  averager,  the  output  of  which  is  plotted  on  one 
channel  of  a  two  channel  chart  recorder .  Figure  3  includes  representative 
Stark  spectra  for  a  range  of  distances  from  the  cathode.  The  splittings  of 
adjacent  components  or  the  width  of  the  entire  Stark  manifold  are  compared  to 
the  Stark  map  of  Fig.  2a  to  determine  the  local  electric  field.  A  relative 
frequency  calibration  is  obtained  by  use  of  an  etalon  with  an  accurately 
determined  spacer.  Interference  fringes  are  generated  and  recorded 
simultaneously  with  the  spectrum. 

The  relative  linestrength  of  Stark  components  or  groups  of  components  is 
also  useful  in  field  measurements . 1 ^  This  approach  works  best  at  low  fields 
where  the  1 inestrengths  are  strongly  dependent  on  field  as  indicated  in 
Fig.  2b.  The  relative  1 inestrengths  are  equal  at  high  field.  In  order  to 
determine  an  electric  field  from  relative  1 inestrengths  one  must  (1)  avoid 
saturation  of  experimental  spectra  and  (2)  establish  that  the  optogalvanic 
effect  is  equally  efficient  on  all  components.  This  second  condition  can  be 
argued  on  fairly  general  grounds.  Collisions  with  ground  state  atoms  rapidly 
mix  the  populations  of  the  zero  field  n  —  11  levels.  The  cross  section  for  8. 
mixing  of  n  «  11  Na  atoms  by  collision  with  ground  state  He  atoms  is 
1  .8rjx10^A‘^.^^  It  is  reasonable  to  assume  that  an  n  =  11  He  atom  behaves 


similarly  to  an  n  =  11  Na  atom,  thus  the  l  mixing  rate  of  approximately 


iixIO^sec  1  at  3.5  Torr  is  expected.  This  mixing  rate  completely  overwhelms 

1  1  _  1 

radiative  decay  from  the  IIP  level  at  a  vacuum  rate  of  1.4x10  sec 

Radiation  trapping  further  strengthens  the  argument  by  lowering  the  effective 

radiative  decay  rate.  A  comparison  of  relative  linestrengths  must  be  based  on 

integrated  linestrengths  because  the  spectral  widths  of  the  Stark  components 

vary  across  the  manifold.  The  outer  components  are  broadened  by  the  local 

electric  field  gradient.  The  laser  probes  some  finite  volume  containing  a 

range  of  field  magnitudes.  The  spectral  position  of  the  outer  components  is 

most  dependent  on  the  field,  and  hence  their  width  is  proportional  to  the  field 

gradient  and  the  volume  probed  by  the  laser. 

The  results  of  field  measurements  as  a  function  of  position  made  at  5 
different  current  densities  in  3-5  Torr  of  He  are  plotted  in  Fig.  4.  These 
measurements  are  based  on  the  splittings  of  adjacent  components  and  the  width 
of  the  Stark  manifold.  The  solid  lines  are  linear  least  square  fits  to  the 
data.  The  linear  decrease  in  field  magnitude  with  distance  from  the  cathode  is 
in  agreement  with  the  classic  picture  of  the  cathode  fall  region.  Close 
inspection  reveals  only  a  very  slight  negative  curvature.  This  linear  behavior 
persists  to  quite  small  fields  near  the  cathode  fall-negative  glow  boundary 
where  the  field  extrapolates  to  zero.  Direct  measurements  to  rather  low  fields 

pq 

have  been  reported  by  Ganguly  et  al .  The  position  of  the  boundary  between  the 
cathode  fall  and  negative  glow  regions  can  also  be  confirmed  by  a  qualitative 
and  quantitative  change  in  optogalvanic  effects.  The  excellent  agreement 
between  the  voltage  measured  using  a  digital  voltmeter  and  the  voltage 
determined  by  integrating  the  electric  field  across  the  cathode  fall  region 
Indicates  that  the  field  measurements  are  on  average  accurate  to  1$. 
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It  is  essential  to  determine  both  the  electric  field  and  the  E/N 
throughout  the  cathode  fall  region.  We  found  that  significant  gas  heating  and 
a  corresponding  gas  density  reduction  occurs  in  the  abnormal  cathode  fall 
region,  even  with  the  water-cooled  electrodes  in  our  experiment.  Ions  carry 
most  of  the  current  in  the  cathode  fall  region.  The  motion  of  atomic  or 
molecular  ions  in  their  parent  gas  (i.e.  He+  in  He  or  Np  in  Np)  is  limited  by 
symmetric  charge  exchange 

He+(fast)  +  He(slow)  •»  He(fast)  +  He+(slow).  (H) 

This  reaction  has  a  large  resonant  cross  section  which  is  only  weakly  dependent 
on  energy.  It  results  in  a  short  equilibration  distance  for  ions  as  discussed 
in  the  next  section,  and  a  very  efficient  conversion  of  electrical  energy  to 
heavy  particle  translational  motion  and  to  heat.  Examination  of  the  typical 
Nd()  products  (gas  density  times  cathode  fall  thickness)  suggests  that  many  of 
the  energetic  neutrals  scatter  before  reaching  the  cathode. 

To  complement  the  field  measurements  in  Fig.  the  gas  temperature  is 
measured  in  the  cathode  fall  region  for  each  of  the  five  current  densities  of 
the  study.  The  temperature  is  found  to  be  constant  within  an  error  bar 
throughout  the  cathode  fall  region. 

Doppler  profiles  are  obtained  by  scanning  a  non-saturating  pulsed  dye 
Laser  through  the  transition  of  interest  and  detecting  the  resulting 
optogalvanic  signal  with  a  boxcar  averager.  The  beam  was  focused  to  a  width  of 
0.06  cm  within  the  cathode  fall  region.  The  transition  is  the  ? ' S  to  3^P 
transition  of  He  at  601.6  nm.  The  dye  laser  bandwidth  is  reduced  to  jOO  Mil/,  by 
introducing  an  etalon.  Other  contributions  to  the  observed  lineshape  include  a 
natural  width  of  92  MHz^®,  pressure  broadening  at  3.3  Torr  of  1H6  MHz^°,  and 
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Stark  broadening  on  the  order  of  100  MHz.  These  contributions  are  small 
compared  to  the  Doppler  width  of  3.66  GHz  at  293  K. 

The  temperature  was  determined  by  measuring  the  full  width  at  half  maximum 
intensity  of  the  line  profile.  Since  the  temperature  increase  is  proportional 
to  the  power  dissipated  in  the  discharge  and  the  Doppler  width  is  proportional 
to  the  square  root  of  the  temperature,  the  square  of  the  measured  width  was 
plotted  as  a  function  of  discharge  power.  If  the  broadening  mechanisms,  apart 
from  Doppler  broadening,  are  negligible  the  resulting  extrapolation  should 
indicate  a  zero  power  Doppler  width  which  corresponds  to  ambient  gas 
temperature.  The  difference  between  the  width  calculated  at  ambient 
temperature  and  that  indicated  by  the  zero-power  extrapolation  gives  the 
residual  width,  which  was  then  subtracted  from  the  measured  widths.  This 
correction  is  less  than  the  size  of  the  error  bar.  The  temperature  at  each 
current  density  was  then  determined  from  the  width. 

Figure  5  is  the  resulting  experimentally  derived  gas  density  as  a  function 
of  discharge  power  measured  at  a  constant  pressure  of  3. 50  Torr.  At  the 
highest  power  studied,  the  density  is  30%  lower  than  what  would  be  expected 
assuming  ambient  temperature.  Each  data  point  is  determined  from  10  Doppler 
width  measurements.  These  density  measurements  are  combined  with  the  field 
measurements  of  Fig.  4  to  derive  important  parameters  of  the  cathode  fall 
region  in  the  next  section  of  this  paper. 

IV.  Current  Balance  at  the  Cathode  Surface 
The  current  balance  at  the  surface  or  the  cathode  is  the  ratio  of  the  ion 
current  to  electron  current,  which  is  determined  from  measurements  of  the 
electric  field  and  gas  density  in  the  cathode  fall  region  and  the  total 
discharge  current  JD.  The  ion  current  density  at  the  cathode,  J°,  is  the 
product  of  the  ion  charge  density,  p+,  and  the  average  ion  velocity  at  the 


cathode  surface,  <vz>°.  Poisson's  equation  determines  the  space  charge  density 
in  the  cathode  fall  region  from  the  spatial  gradient  of  the  electric  field. 

High  fields  in  the  cathode  fall  region  result  in  very  high  electron  velocities 
and  a  negligible  electron  density;  the  space  charge  is  almost  entirely  due  to 
the  ions. 

Ion  motion  in  the  cathode  fall  region  is  limited  by  the  symmetric  charge 
exchange  reaction  of  Eq.  4 .  The  cold  gas  or  high  field  approximation  is 
applicable  throughout  most  of  the  cathode  fall  region.  In  this  approximation 
the  ion  is  assumed  to  start  from  rest  after  each  charge  exchange.  It  is  thus 
reasonable  to  argue  that  the  equilibration  distance  of  the  ions  should  be 
comparable  to  the  mean-free-path  for  symmetric  charge  exchange.  Each  charge 
exchange  collision  wipes  out  all  memory  of  the  ion's  history.  Although  the 
above  argument  is  appealing,  two  aspects  of  the  cathode  fall  region  complicate 
the  picture.  There  is  a  significant  field  gradient,  and  new  ions  are  produced 
throughout  the  cathode  fall  region  by  electron  impact  ionization. 

These  complications  have  both  been  included  in  an  analytic  treatment  of 
ion  transport  in  the  cathode  fall  region.^1  These  solutions  to  the  Boltzmann 
equation  for  ions  will  be  summarized  here.  It  is  convenient  to  reverse 
coordinate  systems  when  discussing  ions.  In  most  of  this  paper,  parameters  are 
measured  and  plotted  as  a  function  of  distance  from  the  cathode.  For  this 
discussion  of  ion  transport  we  choose  an  origin  at  the  boundary  between  the 
cathode  fall  and  negative  glow  regions,  and  we  let  the  positive  /.  axis  point 
toward  the  cathode.  This  temporary  change  in  coordinate  systems  makes  the 
field,  ion  velocity,  and  ion  flux  density  positive.  The  time  independent 
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where  e  is  a  unit  charge,  m+  is  the  ion  mass,  a  is  the  charge  exchange  cross 
section,  and  other  symbols  have  their  usual  meanings.  The  charge  exchange 
cross  section,  although  weakly  energy  dependent,  is  here  approximated  as 
independent  of  energy.  The  first  and  second  terms  on  the  right  hand  side  of 
Eq.  5  are  Wannier's  elegent  expression  for  the  collision  term  due  to  symmetric 
charge  exchange . The  third  term  is  a  source  term  due  to  electron  impact 
ionization.  This  source  term  and  the  electric  field  should  come  from  a 
self-consistent  description  of  the  cathode  fall  region  which  is  beyond  the 
scope  of  this  work.  We  therefore  use  a  field  which  increases  linearly  with 
distance  from  the  cathode  fall-negative  glow  boundary  and  replace  the  ion 
acceleration  eE/m+  by  kz.  The  electron  impact  ionization  source  term,  P(z),  is 
known  to  peak  near  the  cathode  fall-negative  glow  boundary  from  Monte  Carlo 
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simulations  by  Beouf  and  Marode,  and  others.  Section  VII  of  this  paper 
describes  a  determination  of  P(z)  from  Monte  Carlo  simulation  of  electron 
avalanches  in  the  cathode  region.  A  lower  limit  for  the  equilibration  distance 
is  calculated  by  assuming  that  P(z)  describes  a  plane  source  at  the  cathode 
fall-negative  glow  boundary.  An  upper  limit  for  the  equilibration  distance  is 
calculated  by  assuming  P(z)  describes  a  uniform  source  for  positive  z.  Both 
limits  are  derived  from  the  solution  to  the  problem  where  P(z)  is  a  plane 
source  jS(z-z0)  and  where  zQ>0.  The  distribution  function  at  any  positive  z 
for  this  source  is 

fG  =  j5(vx)6(vy)  s(vz)  |exp[-oM(z-z0)]  2  S(vz?  "  kz2  +  kz£)  + 


(6) 
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where  s(vz)  represents  a  step  function  and  other  symbols  have  their  usual 
meanings.  The  integral  ///  fGv  dvz  dvx  dvy  is  the  ion  flux  density,  jz,  for 
all  z  greater  than  zQ.  The  flux  density  must  be  independent  of  position  by 
conservation  of  particles.  The  integral  ///  fG  dvz  dvx  dvy  is  the  ion  density 
which  is  dependent  on  position  because  the  ions  are  accelerating.  The 
distribution  function  of  Eq.  6  for  large  z  approaches  Wannier’s  equilibrium 
distribution  function^2 

feq  =  j6(vx)6(vy)[oN/(kz)]s(vz)exp[-vz2oN/(2kz)].  (7) 

The  lower  limit  for  the  equilibration  distance  of  ions  in  the  cathode  fall 
region  is  determined  by  letting  zQ  approach  zero.  This  limit  corresponds  to  a 
plane  source  of  ionization  at  the  cathode  fall-negative  glow  boundary.  The 
average  velocity  in  the  z  direction  for  the  case  where  zQ  approaches  zero  is 

/^Vzdvzdvxdvy 

<  v  >  =  _ _ 

/;/fGdvzdVxdvy"  ' 
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exp(-oNz)/(/T<z)  +  ttoN  exp(-oNz)[IQ(oNz)+L0(oNz)]/(2A) 


(8) 


where  I0(z)  is  a  modified  Bessel  function  of  order  0  and  LQ(z)  is  a  modified 
otruve  function  of  order  0  as  defined  and  tabulated  in  Abramowitz  and 
Stegun . 33 * 31*  average  velocity  for  this  case  reaches  905t  of  the  equilibrium 
drift  velocity  /2kz/(iroN)  in  1.7  mean-free  paths.  If  zQ>>(oN)  1  then  the  field 
changes  only  slightly  in  a  mcan-free-path  and  the  equilibration  distance  will 


bo  that  of  the  constant  field  problem,  0.65/(oN),  or  two-thirds  of  a 
mean  free-path . 

The  flux  density,  average  velocity,  and  ion  density  are  plotted  in  Fig.  6 
for  the  distribution  function  of  Eq.  6  with  zo=?0.  The  flux  density,  plotted  in 
Fig.  6a,  must  be  independent  of  z  because  the  source  is  localized  at  zQ  near  the 
origin  and  the  system  is  in  steady  state.  The  average  velocity,  plotted  in 
Fig.  6b,  is  a  steadily  increasing  function  of  z  because  the  field  increases 
Linearly  with  z.  Figure  6c  which  is  a  plot  of  the  ratio  of  the  average  velocity 
to  the  local  equilibrium  drift  velocity  J2V.Z/  ltton  ; ,  clearly  shows  that  the 
nonhydrodynamic  behavior  is  confined  to  within  a  few  mean-free-paths  of  the 
source  near  the  origin.  The  equilibration  distance  of  1.7  mean-free-paths 
required  for  the  average  velocity  to  reach  90$  of  the  local  equilibrium  drift 
velocity  is  read  directly  from  Fig.  6c.  A  low  average  velocity  near  the  origin 
implies  a  large  ion  density  because  their  product  is  the  constant  flux  density. 
The  ion  density  is  plotted  in  Fig.  6d.  The  large  ion  density  near  the  origin 
is  due  to  inertia.  A  large  fractional  change  in  ion  density  per  mean-free-path 
is  a  manifestation  of  a  nonhydrodynamic  region.  It  is  because  of  this  large 
density  gradient  that  a  density  gradient  expansion  of  the  continuity  equation 
is  not  applicable  near  the  cathode  fall-negative  glow  boundary.  The  concept  of 
diffusion  must  also  be  abandoned  with  the  density  gradient  expansion.  The 
density  gradient  expansion  and  the  concept  of  diffusion  are  unnecessary  because 
Eq.  6  is  an  exact  solution  to  the  Boltzmann  equation.  These  issues  come  up 
again  later  in  this  paper  when  we  argue  that  there  is  no  "diffusion  like"  ion 
current  across  the  cathode  fall-negative  glow  boundary. 

Equation  6  for  fQ  is  very  slightly  different  from  f^  as  defined  in  Eq.  7 
of  Ref. 31.  The  changes  make  fQ  defined  (zero)  for  positive  z  less  than  zQ. 

Thus  fG  corresponds  more  closely  to  Arfken's  definition  of  a  one-dimensional 
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Green's  function. ^  The  distri bution  function  at  any  positive  z  for  a 
linearly  increasing  field  with  arbitrary  source  function  Is 

f  =  6 (vx)6 (Vy )s(vz)  /  P(zQ) jexp[-oN(z— zQ)]  2 

6 (vz2-kz2+kz02)  +  (oN/Vk)[1-s(vz2-kz2  +  kzQ2)] 

exp{ oNz[/l  -v^/  (kz2 )-1  ]  }  //kz2-vz2 }dzQ.  (9) 


This  integral  determines  a  realistic  ion  distribution  function  If  a  source 
function  from  a  Monte  Carlo  simulation  is  used  as  P(z).  The  objective  here  is 
to  prove  that  the  average  ion  velocity  is  the  equilibrium  drift  velocity  in 
most  of  the  cathode  fall  region.  The  true  source  function  as  determined  in 
Monte  Carlo  simulations  peaks  near  the  cathode  fall-negative  glow  boundary,  and 
thus  is  a  decreasing  function  of  z  throughout  the  cathode  fall  region.  The 
upper  limit  on  the  equilibration  distance  is  determined  by  using  a  uniform 
source  for  all  positive  z.  The  average  ion  velocity  for  a  linearly  increasing 
field  with  a  uniform  source  of  ions  for  all  positive  z  is 


<vz>  - 


Iff  fvzdvzdvxdvy 
77TT  dvzdvxdvy 


exp(-oNz)ir{  I0(oNz)+L0(oNz)  +  oNz[I1  (oNz)+L1  (oNz)  +  2/ir]  }/(2/k) 


where  I^z)  and  L,(z)  are  modified  Bessel  and  modified  Struve  functions  of 
order  one.JJ,J  The  distance  required  for  the  average  velocity  to  reach  90%  of 
the  equilibrium  drift  velocity  /2kz/(iroN)  is  5.7  mean-free-paths. 
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The  ratio  of  the  thickness  of  the  cathode  fall  region  dc,  to  the 
mean-free-path  for  symmetric  charge  exchange  (oN)  1  is  typically  50  to  100  for 
a  rare- gas  cathode  fall.  A  change  in  discharge  current  or  pressure  changes  Ndc 
only  slightly,  thus  the  observation 

oNdc  >>  5.7  (11) 

justifies  using  an  equilibrium  ion  drift  velocity  in  nearly  all  of  the  cathode 
fall  region.  The  approximation  fails  only  within  the  first  6  mean-free-paths 
from  the  cathode  fall-negative  glow  boundary.  The  preceeding  proof  that  the 
ions  are  in  hydrodynamic  equilibrium  enables  us  to  complete  the  determination 
of  the  current  balance  at  the  surface  of  the  cathode. 

The  drift  velocity  at  the  cathode  surface  for  low  currents  is  determined 
from  Helm's  precise  mobility  data.^  At  the  high  current  densities  the  drift 
velocity  is  calculated  using  the  equilibrium  expression 

<vz>°  =  [2eE°/(m+noN)]1/2,  (12) 


where  E°  is  the  field  at  the  cathode  surface.  The  symmetric  charge  exchange 
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cross  section,  o,  is  taken  from  the  calculations  of  Sinha,  Lin,  and  Bardsley. 
These  calculations  agree  with  Helm's  experimentally  derived  cross  sections. 

The  slight  energy  dependence  of  o  does  not  weaken  the  proof  of  equilibration, 
but  it  is  important  to  use  a  value  of  o  that  corresponds  to  the  average  ion 
energy  at  the  cathode  surface.  The  average  ion  energy  at  the  cathode  from  the 
equilibrium  distribution  function  of  Eq.  7  is 


<m+vz2/2>° 


eE°/(2No)  . 


(13) 


A  self-consistent  average  ion  energy,  and  symmetric  charge  exchange  cross 
section  is  conveniently  determined  from  a  graphical  analysis  using  a  plot  of  o 
verses  ion  energy.  The  intersection  of  the  o  verses  ion  energy  plot  with  a 
hyperbola  defined  by  Eq.  13  and  an  empirical  E°/N  determines  the 
self-consistent  average  ion  energy  and  symmetric  charge  exchange  cross  section. 
The  uncertainty  in  <vz>°  is  in  all  cases. 

Table  I  summarizes  the  analysis  of  the  electric  field  and  gas  density 
measurements  to  determine  the  current  balance  at  the  cathode  surface.  The 
excellent  agreement  between  Vvm,  a  voltage  determined  with  a  digital  voltmeter, 
and  vef’  a  voltage  determined  from  integrating  the  field  data,  indicates  that 
the  field  measurements  are  on  the  average  accurate  to  1)f.  The  uncertainty  in 
the  ion  charge  density,  p+,  is  comparable  to  the  uncertainty  in  the  field 
measurements.  The  electron  current  density  at  the  cathode  surface  is  JD-J° 
where  JD  is  the  discharge  current  density.  The  electron  current  density  has, 
of  course,  a  larger  fractional  uncertainty  than  the  ion  current  density.  The 
dimensionless  ratio  J°/(JD~J?)  is  the  current  balance  or  ratio  of  the  ion 
current  to  electron  current  at  the  cathode  surface.  This  ratio  is  of 
particular  interest  because  it  is  the  number  of  ions  produced  in  the  cathode 
fall  region  per  (net)  electron  emitted  from  the  cathode. 

V.  Metastable  Density  Measurements 

The  spatial  dependence  of  2^S  and  2^S  He  metastable  densities  are  mapped 
along  the  axis  of  the  discharge.  The  experiment  has  two  parts:  1) 
laser-induced  fluorescence  is  used  to  create  spatial  maps  of  relative 
metastable  densities  and  2)  absorption  measurements  are  used  to  put  an  absolute 
scale  on  the  density  maps.  The  absolute  density  of  metastables  is  modelled  in 
the  next  section  to  derive  metastable  production  rates  for  comparison  to  Monte 
Carlo  simulations.  From  the  spatial  distribution  of  metastables,  and 
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specifically  from  the  suppression  of  2^S  metastables  in  the  negative  glow, 
valuable  information  is  gained  about  the  low  energy  electrons  in  the  negative 
glow . 

Relative  metastable  density  maps  are  produced  for  each  of  the  5  current 
densities  under  study.  The  laser-induced  fluorescence  measurements  are  made 
using  a  pulsed  N2  laser  pumped  dye  laser  tuned  to  either  the  501  .6  nm 
transition  for  mapping  2^  density  or  the  388.9  nm  transition  for  mapping  the 
2^s  density.  The  beam  passes  through  a  250  pm  slit  which  was  imaged  1:1  into 
the  center  of  the  discharge  using  a  1 3  cm  focal  length  cylindrical  lens.  The 
laser  polarization  is  along  the  axis  of  the  discharge.  Fluorescence  is 
collected  in  a  direction  mutually  perpendicular  to  the  laser  axis  and  the  axis 
of  the  discharge.  A  lens  and  slit  combination  is  used  to  isolate  the  central 
portion  of  the  laser-discharge  interaction  region.  Thi3  eliminates  possible 
distortion  of  the  spatial  map  due  to  meta3table  populations  in  the  fringes. 

The  lens  is  masked  to  cut  down  on  scattered  fluorescence  which  can  also  distort 
the  map. 

Excitation  transfer  among  the  n  -  3  levels  due  to  collisions  with  ground 
state  atoms  allows  the  monitoring  of  fluorescence  at  a  wavelength  different 
from  that  of  the  laser. 3®  Interference  filters  are  used  to  isolate  667.8  nm 
( 3 1 D  -  2^ P )  radiation  when  the  laser  is  tuned  to  501.6  nm  and  706.5  nm 
(3^S  -  2^p)  radiation  when  the  laser  is  tuned  to  388.9  nm.  In  this  way 
scattered  laser  light  is  eliminated  from  the  signal.  The  fluorescence  is 
detected  with  a  photomultiplier.  The  entire  fluorescence  collection  system  is 
mounted  on  the  translation  stage  so  it  remains  fixed  with  respect  to  the 
discharge.  The  signal  is  digitized  and  integrated  using  a  Transiac  2001S 
transient  recorder  with  a  Ml 00  averaging  memory  interfaced  with  an  IBM  PC. 
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The  laser  absorption  measurements  for  determining  an  absolute  scale  for 
the  metastable  densities  are  made  using  the  same  transitions  as  the 
laser-induced  fluorescence  studies.  An  etalon  is  introduced  to  reduce  the 
laser  bandwidth  to  500  MHz.  The  beam  passes  through  a  500  pm  slit  which  is 
imaged  1:1  into  the  discharge.  Absorption  measurements  are  made  at  the  peak  of 
the  spatial  density  at  two  current  densities.  The  intensity  of  the  laser  is 
substantially  reduced  to  avoid  saturation  of  the  transition.  The  transmitted 
laser  light  is  monitored  with  a  photodiode  detector.  The  signal  from  the 
detector  is  processed  with  a  boxcar  averager  and  plotted  on  a  strip-chart 
recorder  as  the  laser  is  scanned  through  the  transition.  Metastable  densities 
are  determined  from  the  integral  over  frequency  of  the  natural  log  of  the 
transmittance.  Each  density  is  an  average  of  10  measurements. 

Figures  7  and  8  show  the  results  of  the  metastable  density  measurements. 
The  density  maps  have  an  uncertainty  of  -20%  in  the  center  of  the  discharge 
rising  to  -30$  near  the  electrodes. 

VI.  Metastable  Transport  and  Kinetics 

The  transport  and  kinetics  of  2'S  and  2^s  metastables  in  the  discharge  is 
modelled  using  a  pair  of  balance  equations  which  are  coupled  differential 
equations: 
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Ms  t  =  metnstable  density 

N  =  density  of  ground  state  atoms 
ng  =  density  of  low  energy  electrons 
D  =  diffusion  coefficient 
P  =  production  rate  per  unit  volume 

Y  =  rate  constant  for  singlet  metastable  destruction  due  to  collisions  witl 
ground  state  atoms 

6  =  rate  constant  for  destruction  of  metastables  due  to  metastable-metasta| 

i 

collisions 


<  =  effective  rate  constant  for  destruction  of  singlet  metastables 
due  to  low  energy  electron  collisions 
s  and  t  subscripts  indicate  singlet  and  triplet  meta3tables  respectively. 
The  first  three  terms  in  each  equation  are  lo3S  terms  arising  from  diffusion 
and  collisions  between  metastables.  There  is  an  additional  loss  term  in  the 
singlet  equation  due  to  collisions  with  ground  state  He  atoms.  The  marked 
suppression  of  singlet  metastables  in  the  negative  glow  (refer  to  Fig.  7)  is 
due  to  low  energy  electron  collisions,  i.e.  the  reactions: 


He(21S)  +  e  He(23S)  +  e  +  0.79  eV 
-  He(21P)  +  e"  -  0.60  eV. 


(16a) 

(16b) 


These  reactions  both  have  large  cross  sections.  Exothermic  reaction  16a  is  a 
spin  exchange,  s-wave  process  with  a  large  cross  section  for  low  energy 
electrons.  Endothermic  reaction  16b  corresponds  to  a  very  strong  optically 
allowed  transition;  it  has  a  large  cross  section  peaking  at  electron  energies 
roughly  three  times  threshold  energy.  Reaction  16a  dominates  the  2^ 
quenching,  because  1 6b  must  compete  with  the  2^P  -*  2^  radiative  decay.  The 


metastable  spin  conversion  reaction  I6ais  included  as  the  fifth  loss  term  in 
the  singlet  equation  and  a  corresponding  gain  in  the  triplet  equation. 

Low  energy  electrons  are  confined  to  the  negative  glow  region  due  to  the 
high  fields  in  the  cathode  fall  region.  The  spatial  dependence  of  the  low 
energy  electron  density,  ng ,  is  therefore  that  of  a  step  function;  it  is  zero 
in  the  cathode  fall  region  rising  sharply  at  the  negative  glow  boundary  to  some 
constant  value.  The  spatial  dependence  of  the  production  terms  is  assumed  to 
be  a  fundamental  diffusion  mode  except  at  the  highest  and  lowest  currents  where 
the  fundamental  mode  is  skewed  slightly  to  maintain  the  peak  at  the  cathode 
fall-negative  glow  boundary.  The  assumption  that  the  production  peaks  at  the 
boundary  i3  borne  out  by  Monte  Carlo  simulations  discussed  later  in  this  paper. 
Values  for  DN,  Y  and  8  are  taken  from  Phelps. ^  The  temperature  dependence  of  D 
is  taken  from  Buckingham  and  Dalgarno,140  that  of  Y  is  taken  from  Allison, 

Browne  and  Dalgarno. 

The  key  quantities  which  are  determined  in  this  modelling  are  the 
integrated  production  rates,  / Pgdz  and  /Ptdz,  and  the  product  of  the  rate 
constant  for  metastable  spin  conversion  with  the  low  energy  electron  density, 
<ng.  Values  for  these  key  quantities  are  assumed,  the  coupled  equations  are 
solved  for  Ms(z)  and  Mt(z),  and  a  comparison  is  made  to  the  empirical 
metastable  density  maps.  Diffusion  modes  up  to  tenth  order  are  included  in  the 
solutions  because  of  the  spatial  assymmetry  introduced  by  metastable  spin 
conversion.  After  the  comparison  with  the  empirical  metastable  maps,  the 
values  for  /Pgdz,  /Ptdz,  and  *:ng  are  adjusted  and  the  processes  iterated  until 
a  good  fit  to  the  empirical  metastable  maps  is  achieved.  The  value  for  <ne 
determined  in  this  fashion  should  be  considered  an  effective  rate  of  metastable 
spin  conversion  because  it  must  be  modified  to  include  the  reverse  reaction. 

The  reverse  reaction  is  not  negligible  if  the  relative  singlet  verses  triplet 
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metastable  density  approaches  thermal  equilibrium  with  the  low  energy  electron 
gas  in  the  negative  glow. 

Figures  9  and  10  show  both  the  measured  and  calculated  2^  and  2% 
densities  respectively  for  one  of  the  current  densities.  Table  II  summarizes 
for  all  current  densities  the  parameters  used  in  the  model  such  as  Ds,  Dt,Y, 
and  B  and  the  key  quantities  /Pgdz,  /Ptdz,  and  <ne.  The  dimensionless  ratios 
e/P  dz/ ( JD~J°)  and  e/Ptdz/ ( Jq-J° ) ,  where  Jp-J°  is  the  empirical  electron 
current  density  at  the  cathode,  are  of  particular  interest.  These 
dimensionless  ratios  are  the  average  number  of  singlet  and  triplet  metastables 
produced  per  (net)  electron  emitted  from  the  cathode. 

VII.  Monte  Carlo  Simulations 

Monte  Carlo  simulations  are  now  widely  used  to  study  the  behavior  of 
charged  particles  in  weakly  ionized  gases.  These  simulations  are  not  as 
elegant  as  a  numerical  solution  of  the  Boltzmann  equation,  and  they  do  require 
considerably  more  computer  time  than  numerical  solutions.  The  computer  time 
issue  is  steadily  decreasing  in  significance;  all  of  the  simulations  described 
in  this  work  were  executed  on  a  personal  computer.  Monte  Carlo  simulations  are 
very  easy  to  code,  and  they  can  be  made  quite  realistic  if  sufficient  high 
quality  data  is  available  on  cross  sections  including  angular  distributions  for 
elastic  scattering,  for  excitation,  and  for  ionization.  It  is  straightforward 
to  include  almost  any  amount  of  detail  in  a  simulation  which  is  justified  by 
the  data  base  of  cross  sections.  Any  well  defined  physical  quantity  can  be 
determined  in  a  Monte  Carlo  simulation.  Monte  Carlo  simulations  are  often  used 
to  test  numerical  solutions  of  the  Boltzmann  equation  which  usually  require 
some  approximations.  The  simulations  are  also  very  useful  for  comparison  to 
experiments . 


Monte  Carlo  simulations  are  used  in  this  work  to  study  electron  avalanches 
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in  the  cathode  fall  and  negative  glow  regions.  There  is  no  need  for  Monte 


Carlo  simulations  of  ion  transport  in  the  cathode  fall  region  because  we  have 


the  exact  analytic  solutions  to  the  Boltzmann  equation  for  ions  as  presented  in 


Section  IV.  The  Monte  Carlo  simulations  include  the  empirical  field 


distributions  and  gas  densities  for  the  cathode  fall  region  as  given  in  Section 


III.  A  uniform  field  of  1.0  V/cm  or  10.0  V/cm  is  assumed  for  the  negative  glow 


region.  We  suspect  that  there  is  actually  a  reversal  of  the  field  direction  at 


the  boundary  between  the  cathode  fall  and  negative  glow  regions.  The  reversal 


would  produce  a  potential  energy  well  which  holds  the  low  energy  electron  gas 


of  the  negative  glow.  The  reversal  would  also  result  in  a  back  diffusion  (or 


drift)  toward  the  anode  of  ions  produced  in  the  negative  glow.  Unfortunately 


it  is  difficult  to  include  a  field  reversal  in  the  Monte  Carlo  simulations. 


However,  the  actual  magnitude  and  direction  of  the  weak  field  in  the  negative 


glow  region  has  very  little  effect  on  excitation  and  ionization  in  the  negative 


glow  which  is  produced  by  high  energy  beam  electrons  from  the  cathode  fall 


region . 


The  Monte  Carlo  simulations  are  three  dimensional  in  the  sense  that  they 


include  angular  scattering,  but  they  do  not  include  edge  effects  due  to 


fringing  of  the  discharge.  One  dimensional  infinite  plane  parallel  geometry  is 


assumed.  The  most  convenient  coordinates  are  T,  y,  and  z  where  T  is  the 


electron's  kinetic  energy,  y  is  the  cosine  of  the  angle  between  the  electron's 


velocity  and  acceleration,  and  z  is  the  distance  from  the  cathode.  The 


acceleration  5  of  magnitude  |eE/me|  is  in  the  positive  z  direction.  We  ask  the 


reader  to  note  that  we  have  returned  to  our  original  coordinate  system,  which 


was  abandoned  only  during  the  discussion  of  ion  transport  in  Section  IV.  The 
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quantity  s  is  a  distance  along  an  electron's  trajectory  between  collisions,  and 
W  is  electron's  potential  energy. 

Two  quantities  are  conserved  while  the  electron  is  in  free  flight  between 
collisions,  the  transverse  kinetic  energy  T*  where 


T*  =  T(1-/) 


and  the  total  energy  T  +  W(z).  An  electron  is  launched  after  a  collision  with 
initial  conditions  T1 ,  and  z1 .  Once  a  single  chosen  coordinate  is 
determined  immediately  before  the  next  collision,  such  as  Tc,  uc,  or  zc,  the 
other  two  coordinates  are  determined  from  conservation  equations.  The  change 
in  the  chosen  coordinate  is  determined  from  a  pseudo  random  number  in  a  fashion 
that  reflects  the  electron's  scattering  probability  along  its  trajectory.  By 
tracking  large  numbers  of  electrons  one  can  determine  the  same  information 
embodied  in  a  solution  to  the  Boltzmann  equation:  a  phase  distribution  function 
and  all  moments  of  it.  Tne  enormous  advantage  of  a  null  collision  Monte  Carlo 
technique  is  that  it  makes  it  possible  to  relate  the  change  in  the  chosen 
coordinate  between  collisions  to  the  pseudo  random  number  without  performing  a 
numerical  integration  of  scattering  probability  along  the  trajectory. 

Numerical  integration  is  a  slow  process  and  is  to  be  avoided  if  possible. 

The  null  collision  technique  used  in  this  investigation  was  developed  by 
Boeuf  and  Marode  for  nonuniform  fields.  1  The  probability  density  function 
q(sc)  is  used  to  describe  the  differential  probability  of  a  collision  oocuring 
in  a  interval  ds  about  sc, 


dp  =  q( s  )ds  -  exp{-/  c  No[T(s)]ds}  No[T(s.)]ds. 
c  o  c 


The  electron  starts  with  s-0  after  the  last  collision  at  coordinates  T1 ,  v1 , 


The  cross  section  a  includes  elastic,  excitation,  and  ionization 


scattering.  It  is  obviously  very  dependent  on  the  electron  kinetic  energy  T. 
The  exponential  in  the  definition  of  q  is  the  probability  of  the  electron 
surviving  without  a  collision  to  sc,  and  the  remaining  factor  is  the 
probability  of  scattering  at  sc.  The  cumulative  distribution  function 


P(sc) 


;s° 

o 


s  s 

dp  =  /  c  q(s)ds  =  1-exp{-/Qc  No[T(s)]ds} 


(19) 


is  a  probability  of  a  collision  occuring  before  the  electron  reaches  sc-  This 
function  should  be  set  equal  to  a  pseudo  random  number  r.j  on  the  interval  from 
0  to  1 .0  and  solved  for  sc  to  determine  Tc  or  uc  or  zc  in  the  Monte  Carlo 
simulation , 

g 

/  0  No[T ( s) ]ds  -  -lnd-r,)  =  R.  (20) 

o 

Following  Boeuf  and  Marode  we  note 

ds  =  dz/u  (21 ) 


and 


adz  =  dT 


(22) 


thus 


;Tc  No(T) 
T.  a  u 


/c  No(T)  /T 
T1  a(T-T«)1/2 


sgn(y)dT 


R. 


(23) 


Note  that  the  acceleration  a,  which  is  a  function  of  position,  must  be 
expressed  as  a  function  of  T  in  order  to  evaluate  this  integral.  Of  course  the 
null  collision  technique  avoids  this  numerical  integration.  There  is  an  upper 
bound,  B^,  of  No(T)/T/a  in  the  ith  spatial  bin  defined  by  the  interval  z^  to 
zi+1.  The  integration  of  Eq.  23  is  straightforward  if  a  null  scattering 
process  is  included,  with  a  cross  section  chosen  such  that 

N[o(T)  ♦  °null (T ) ]  /T/a  =  Bi  (2*0 

throughout  the  ith  bin  from  z^  to  z<  +  1.  Null  collisions,  when  they  occur,  do 
not  change  the  electron's  velocity.  The  remaining  complications  in  the 
integral  of  Eq.  23  are  due  to  sgn(p)  and  to  the  steps  in  the  integrand  caused 
by  Bi_1  1  . . . .  These  complications  are  handled  by  appropriate  branches 

and  loops  in  the  program.  First  we  consider  a  non-negative  •  The  equation 

2Bi(/TrT*  -  /T1 -T* )  =  R  (25) 

is  solved,  and  the  solution  T  is  compared  to  T^+1,  the  kinetic  energy  at  the 
boundary  between  the  i '  th  and  i'th+1  bin  at  z^+i*  If  T  is  less  than  Tj*+i»  then 
T=TC  and  a  collision  occurs  in  the  i'th  bin.  Otherwise  the  equation 

2Bi+1 (/T-T*  -  /Tf+1-T»)=R-2Bi(/T5+1rT*  -  /T 1 ~T*)  (26) 

is  solved,  and  the  solution  T  is  compared  to  T^+2-  This  procedure  is  iterated 
until  the  bin  and  the  kinetic  energy  Tc  immediately  before  the  collision  is 
determined.  A  similar  procedure  is  followed  for  negative  wi  with  the 
additional  possibility  of  the  electron  reversing  direction  before  colliding. 


-35- 


F or  exampl o  we  consider  an  electron  with  negative  p ,  but  without  sufficient 
longitudinal  kinetic  energy  T.-T*  to  reach  zi .  The  equation 

-2Bi(/T-T*  -  /T ■,  -T* )  =  R  (27) 

can  be  solved  for  T  only  if  R  is  sufficiently  small.  A  large  value  of  R 
implies  that  the  electron  passes  through  the  turning  point  where  T=T*  and  p=0; 
thus  the  equation 

2B.(/T-T*  +  /Tj-T* )  -  R  (28) 

is  solved  for  T. 

The  type  of  collision  is  determined  by  choosing  a  second  pseudo  random 
number  r2  in  the  interval  from  0  to  1.0.  The  0  to  1.0  interval  is  divided  into 
segments,  each  of  length  equal  to  the  probability  of  a  particular  type  of 
collision  at  energy  Tc.  The  cross  section  o  is  a  sum  of  the  elastic  cross 
section  oel ,  the  ionization  cross  section  Cj;on>  and  various  excitation  cross 
sections.  The  first  segment  is  of  length  0nun (TQ )/[onull (Tc)  4  o(Tc)],  the 
second  is  of  length  <Jel(Tc)/[onull(Tc)  +  o(Tc)],  the  third  is  of  length 
°ion^Tc^°null  ^c'*  +  and  so  on  f’or  each  excitation  cross  section.  The 

segment  containing  r2  determines  the  type  of  collision. 

The  electron's  coordinates  T2  an<3  \*2  *mmediately  after  the  collision  must 
be  determined.  The  procedure  is  dependent  on  the  type  of  collision.  Null 
collision,  of  course,  result  in  T2  =  Tc  and  p2  =  uc.  The  determination  of  T2 
and  p2  for  real  collisions  requires  differential  cross  section.  The  sources  of 
the  cross  sections  used  in  the  simulation  are  given  during  the  discussion  of  T2 
and  up* 
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The  total  and  differential  elastic  scattering  cross  section  for  electron 

|l  O 

impact  of  He  calculated  by  LaBahn  and  Callaway  are  used  in  the  simulations. 

The  coordinate  v?  after  an  elastic  collision  is  calculated  from  the  scattering 
angle  x  and  the  azimuthal  scattering  angle  $  using  the  equation 

u2  =  uccosx  +  /I sinx  cos<j>.  (29) 


The  scattering  angles  x  and  4>  are  determined  from  two  more  pseudo  random 
numbers  using  the  equations 


r 


3 


2ir  ;X  dqel 
°el  o  dfl 


sinGdG 


(30) 


and 


r4  =  <J>/ 2ir. 


(31 ) 


Equation  30  can  be  solved  for  x  by  numerical  integration  of  LaBahn  and 

li  p 

Callaway's  differential  cross  section.  Approximations  such  as 


— -r~  =  A[1+g(Tc)sin2(X/2)]_1 


(32) 


are  particularly  convenient  because  x  can  be  evaluated  without  numerical 
integration.  The  function  g(TQ)  is  roughly  quadratic  in  Tc •  The  kinetic 
energy  after  the  collision  is 

T2  =  Tc  Cl-(2rae/m0)(1-cosx)],  (33) 
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whi-r*'  inu  is  the  atom  mass. 


The  analytic  expressions  given  by  Alkhazov  for  excitation  cross  sections 


are  used  in  the  simulations.  J  Excitation  of  n«2  through  n=5  singlet  and 


triplet  s  and  p  levels  is  included  in  the  simulation.  Excitation  of  the  n=*3 


through  n=5  singlet  and  triplet  d  levels  is  also  included.  There  is  very 


little  experimental  or  theoretical  information  on  differential  excitation  cross 


sections  for  electrons  on  He.  Thus  we  use  an  isotropic  scattering 


approximation  for  excitation, 


ii2  =  2r^  -  1  .0 


The  kinetic  energy  immediately  after  a  excitation  collision  is 


T2  =  Tc  ~  “exc 


where  We  is  the  appropriate  excitation  energy.  Recoil  is  neglected  in  this 


approximation . 


A  fraction  of  the  atoms  excited  to  the  3p  levels  and  above  are 


associatively  ionized.  An  accurate  determination  of  the  fraction  requires 


accurate  associative  ionization  cross  sections,  cross  sections  for  excitation 


transfer  due  to  ground  state  atom  collisions,  and  radiative  decay  rates.  On 


the  basis  of  what  information  is  available  we  estimate  that  25 %  of  the  atoms  in 
the  3p  levels  and  above  are  associatively  ionized  at  3.5  Torr,^®*^®  TV,° 


electrons  released  in  associative  ionization  contribute  to  the  avalanche,  they 


are  assumed  to  have  an  initial  kinetic  energy  of  1 .0  eV  and  an  isotropic 


angular  distribution.  The  remaining  75%  of  the  atoms  excited  to  the  3P  and 


higher  triplet  levels,  and  all  of  the  atoms  excited  to  the  3^S  and  2^P  levels 
are  assumed  to  radiate  to  the  2^S  metastable  level.  This  cascade  contribution 
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is  combined  with  direct  excitation  of  the  2^S  level  to  determine  the  total 
production  rate  for  2^S  metastables.  The  total  production  rate  for  2^ 
metastables  includes  direct  excitation  and  a  cascade  contribution  of  19%  of  the 
excitation  to  the  3p  and  higher  singlet  levels.  The  total  production  rate  for 
the  2^?  resonant  level  includes  direct  excitation  and  a  cascade  contribution  of 
all  of  the  3^  excitation  and  56%  of  the  excitation  to  the  3p  and  higher 
singlet  levels.  These  assumptions  are  based  on  a  significant  l  mixing  of  3p 
and  higher  level  populations  due  to  collision  with  ground  state  atoms. 

In  a  vacuum  nearly  all  of  the  atoms  excited  to  the  2^P  resonant  level 
decay  to  the  ground  state  via  VUV  emission.  Radiation  trapping  reduces  the 
effective  decay  rate  to  1  .3x10^sec_1  ,  which  is  comparable  to  the  vacuum  decay 
rate  of  1  .98x1 O^sec-1  for  the  21P->21S  infrared  branch  at  2um.2®  However  the 
21P-»21S  branch  is  also  trapped  in  the  cathode  fall  region  of  our  discharge 
experiments,  thus  we  do  not  include  any  of  the  2^P  excitation  in  calculating  a 
total  production  rate  for  2^  metastables.  This  approximation  is  not 
completely  satisfactory  because  the  2^P  ■*  2^S  branch  is  not  highly  trapped  in 
the  negative  glow  where  the  2^S  density  is  suppressed  by  low  energy  electron 
collision.  The  21?  level  is  feeding  population  to  the  2^  level  through  a 
col  1 isional-radiati ve  coupling,  and  the  2^S  level  is  feeding  population  to  the 
2%  through  a  collisional  coupling.  The  net  result  is  that  the  Monte  Carlo 
simulation  produces  a  metastable  production  rate  which  is  somewhat  lower  than 
the  empirical  rate. 

The  analytic  expression  by  Alkhazov  for  the  total  ionization  cross  section 
is  used  in  the  simulation . ^  Alkhazov  also  provides  an  analytic  expression  for 
the  differential  cross  section  with  respect  to  energy. ^  The  kinetic  energy  of 
the  outgoing  scattered  electrons  is  determined  from  a  pseudo  random  number 
using  the  expression, 
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J_  /2  f!*2S  dT. 


°ion  ° 


dT 


(36) 


The  outgoing  ejected  electron  has  kinetic  energy 


T2  ”  Tc  ‘  T2  wion 


(37) 


where  wion  is  the  ionization  potential.  The  coordinates  u2  and  Ug*  are 
determined  one  of  two  ways. 

The  first  simulations  assumed  isotropic  angular  distributions  for  both 
outgoing  electrons  in  an  ionizing  collision.  Some  of  the  results  based  on  the 
isotropic  angular  distribution  for  ionizing  collisions  were  presented  in  a 

UD 

recent  Phys.  Rev.  Letter.  The  first  simulations  predicted  a  current  balance 
in  good  agreement  with  experiment  in  a  near  normal  cathode  fall  region,  but 
predicted  too  much  ionization  in  a  highly  abnormal  cathode  fall  region.  The 
E/N  reaches  3000Td  in  the  abnormal  cathode  fall  region,  thus  the  inelastic 
processes  are  resulting  in  forward  peaked  angular  scattering  distributions. 
Numerical  solutions  of  the  Boltzmann  equation  for  electrons  have  traditionally 
assumed  isotropic  angular  distributions  for  inelastic  collisions.  Recent  work 
by  Phelps  and  Pitchford  explores  the  effect  of  anisotropic  inelastic  scattering 
at  very  high  E/N  in  N2.^ 

Our  second  simulations  use  the  same  assumptions  used  by  Boeuf  and  Marode: 
(1)  the  incident,  scattered,  and  ejected  electron  velocities  are  coplanar,  and 
{?.)  the  scattered  and  ejected  electron  velocities  are  perpendicular.1^  These 
assumptions  result  in 


cosx  -  /T2/(Tc  -  Wion) 


(38) 
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and 

cosX '  =  /T2*/(Tc-Wion)  .  (39) 

The  azimuthal  scattering  angle  $  is  determined  by  a  pseudo  random  number  as  in 
Eq.  31 .  Equation  29  is  used  to  determine  The  coordinate  calculated 

using 


u2'  =  Pc  cosx'  -  /l-yc2  sinx’  cos<{>.  (40) 

The  justification  for  this  angular  distribution  is  not  strong,  but  it  does 
represent  a  reasonable  attempt  to  introduce  anisotropic  scattering  in  ionizing 
collisions. 

The  details  of  the  electron  interactions  with  the  cathode  and  anode  must 
be  specified.  The  kinetic  energy  of  electrons  emitted  from  the  cathode  is 
assumed  to  be  5.0  eV  and  they  are  assumed  to  have  a  random  angular  distribution 
in  the  forward  direction.  The  use  of  a  more  realistic  energy  distribution  for 
electrons  emitted  at  the  cathode  does  not  affect  ionization  or  excitation  in 
the  Monte  Carlo  simulations.  Some  electrons  are  reflected  back  to  the  cathode 
by  elastic  scattering  from  He  atoms  before  they  suffer  an  inelastic  collision. 
Electrons  do  not  have  sufficient  total  energy  to  reach  the  cathode  after  an 
inelastic  collision.  The  percentage  of  the  gross  electron  emission  which  is 
back  scattered  to  the  cathode  varies  from  19$  at  the  lowest  current  density  (a 
near  normal  cathode  fall)  to  3.6$  at  the  highest  current  density  (a  highly 
abnormal  cathode  fall).  All  Monte  Carlo  results  are  presented  in  terms  of  net 
electron  emission  from  the  cathode.  The  Monte  Carlo  simulation  assumes  no 
reflection  from  the  anode,  each  electron  is  absorbed  immediately  when  it 
reaches  the  anode.  This  approximation  will  cause  the  Monte  Carlo  simulation  to 
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predict  excitation  and  ionization  rates  in  the  negative  glow  which  are  lower 
than  in  experiment.  The  discordance  is  significant  only  in  the  case  of  the 
highly  abnormal  glow  discharge  at  1.50  mA /cm  . 

Figure  11  includes  Monte  Carlo  histograms  giving  the  number  of  ionization 

and  excitation  events  per  (net)  electron  emitted  from  the  cathode  as  a  function 

of  distance  from  the  cathode.  The  results  of  Fig.  11  are  based  on  the 

anisotropic  angular  distribution  for  ionizing  collisions.  These  histograms  are 

for  field  distributions  and  gas  densities  corresponding  to  0.519  mA/cm2.  The 

spatial  dependence  of  the  ionization  and  excitation  rates  are  all  described  by 

a  roughly  symmetric  function  which  peaks  near  the  cathode  fall-negative  glow 

boundary.  The  simulations  show  that  the  "beam"  electrons  in  He  penetrate 

beyond  the  cathode  fall-negative  glow  boundary  by  a  distance  comparable  to  d  , 

c 

the  thickness  of  the  high  field  cathode  fall  region.  The  penetration  distance 
will  scale  as  the  inverse  of  the  gas  density  in  a  fashion  similar  to  dc.  The 
'•  S  metastable  excitation  has  the  most  asymmetric  spatial  dependence;  it  is 
stronger  in  the  negative  glow  than  in  cathode  fall  region.  This  i3  explained 
by  the  sharp  peak  just  above  threshold  in  the  energy  dependence  of  triplet 
excitation  cross  sections.  The  asymmetry  in  the  empirical  maps  of  the  2% 
metastable  density  is  due  in  part  to  the  spatial  dependence  of  triplet 
excitation.  The  spatial  dependences  of  excitation  and  ionization  at  other 
discharge  current  densities  are  similar.  The  results  of  Monte  Carlo 
simulations  for  all  five  current  densities  are  presented  in  Tables  III  and  IV. 
The  Table  III  results  are  based  on  the  isotropic  angular  distribution  for 
ionizing  collisions,  and  those  of  Table  IV  are  based  on  the  anisotropic  angular 
distribution  for  ionizing  collisions. 


The  absolute  accuracy  of  the  Monte  Carlo  simulations  must  be  estimated  in 

order  to  compare  the  simulations  to  experiment.  Statistical  uncertainty  is  not 

large  because  several  thousand  complete  avalanches  were  simulated  for  each  of 

the  five  discharge  currents.  Alkhazovs'  analytic  expressions  for  the  He  cross 
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sections  are  based  on  theoretical  and  experimental  work  by  many  groups.  J  He 
estimates  that  the  uncertainty  is  ±25?  at  low  energy,  and  ±5?  at  high  energy 
where  a  Born  approximation  is  reliable.  We  believe  his  uncertainties  are 
reasonable  and  we  estimate  that  the  Monte  Carlo  results  have  an  uncertainty  of 
±25?. 

VIII.  Comparison  of  Experimental  Results  and  Monte  Carlo  Simulations 

The  first  comparison  of  experimental  results  to  Monte  Carlo  simulations 

involves  the  current  balance  at  the  cathode.  This  dimensionless  quantity  is 

the  ratio  of  ion  to  electron  current  at  the  cathode.  Figure  12  is  a  plot  of 

the  ratio  versus  total  discharge  current  density.  The  experimental 

measurements  are  points  with  uncertainties  which  are  discussed  in  Section  IV. 

The  dotted  line  is  the  same  ratio  determined  from  Monte  Carlo  simulations  using 

the  isotropic  angular  distribution  for  ionizing  collisions.  The  calculations 
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represented  by  the  dotted  line  were  included  in  a  recent  Phys.  Rev.  Letter. 

The  solid  line  is  the  ratio  as  determined  from  Monte  Carlo  simulations  using 
the  anisotropic  angular  distribution  for  ionizing  collisions.  All  Monte  Carlo 
simulations  include  the  same  anisotropic  elastic  scattering.  Clearly  the 
assumption  of  isotropic  Inelastic  scattering  is  beginning  to  fail  in  the  highly 
abnormal  cathode  fall  region  where  E/N  reaches  3000  Td. 

One  of  the  key  issues  in  interpreting  the  Monte  Carlo  results  is  whether 
or  not  one  should  include  any  of  the  ionization  produced  in  the  negative  glow 
when  computing  the  ratio  of  the  ion  to  electron  current  at  the  cathode.  We 
maintained  in  a  recent  Phys.  Rev.  Letter,  and  we  continue  to  maintain  that 


essont. i  al  ly  none  of  the  ionization  produced  in  the  negative  glow  contributes  to 

il 

the  ion  current  at  the  cathode.  The  ratio  of  the  ion  to  electron  current  at 
the  cathode  is  computed  by  summing  ionization  in  spatial  bins  of  a  Monte  Carlo 
histogram  such  as  Fig.  11.  The  sum  includes  spatial  bins  from  the  cathode  to  a 
distance  dc  from  the  cathode  where  the  cathode  fall  -  negative  glow  boundary 
occurs.  The  justification  for  neglecting  any  "diffusion"  of  ions  across  the 
boundary  is  postponed  until  after  a  discussion  of  the  negative  glow. 

The  current  balance  from  the  Monte  Carlo  simulations  includes  direct 
electron  impact  ionization  and  associative  ionization  due  to  collisions  between 
highly  excited  and  ground  state  He  atoms.  Ionization  from  collisions  between 
pairs  of  metastable  atoms  is  not  included,  but  it  is  a  3mall  (£10%) 
contribution  to  the  current  balance  from  the  Monte  Carlo  simulations. 

The  current  balance,  typically  3-3,  is  lower  than  we  first  expected.  It 
is  also  surprisingly  independent  of  total  discharge  current.  The  greater 
electric  field  and  thickness  in  the  cathode  fall  region  at  high  current  are 
offset  by  a  gas  density  reduction  due  to  heating.  The  net  result  is  that  on 
average  an  electron  produces  the  same  amount  of  ionization  in  the  cathode  fall 
region  at  low  and  high  discharge  currents.  It  is  not  yet  possible  to  say 
whether  a  current  balance  of  3*3  is  typical:  for  only  He,  for  inert  gases,  or 
for  many  gases. 

The  second  comparison  involves  the  number  of  metastable  atoms  produced  in 
the  cathode  fall  and  negative  glow  regions  per  (net)  electron  emitted  from  the 
cathode.  The  experimental  measurements  are  plotted  as  points  in  Fig.  13.  The 
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empirical  metastable  production  in  units  of  cm  3ec  ’  is  divided  by  the 
empirical  electron  flux  density  at  the  cathode  (Jp-J°)/e.  The  error  bars 
reflect  onLy  uncertainties  in  the  metastable  density  measurements  a3  discussed 
in  Section  V.  The  solid  line  of  Fig.  13  is  the  result  of  Monte  Carlo 


simulations  using  the  anisotropic  angular  distribution  for  ionizing  collisions. 


The  21 S  production  plotted  in  Fig.  13a  shows  excellent  agreement  between 
experiment  and  the  Monte  Carlo  simulations  except  at  0.190  mA/cm^.  The 
discordance  in  the  near  normal  glow  discharge  at  0.190mA/cm2  is  in  part  due  to 
the  neglect  of  the  2^?  ■*  2^  cascade  contribution  in  Monte  Carlo  determinations 
of  2 's  production.  At  low  current  densities  the  21  P  -*  2'*S  radiation  is  not 
completely  trapped.  The  experimental  2^S  production  plotted  in  Fig.  13b  is 
somewhat  higher  than  that  predicted  by  the  Monte  Carlo  simulations.  The  total 
2  S  and  2JS  metastable  production  plotted  in  Fig.  13c  is  also  somewhat  higher 
than  that  predicted  by  the  Monte  Carlo  simulations.  At  least  two  effects  are 
responsible.  The  experimental  2^S  production  and  the  total  metastable 
production  are  enhanced  by  a  coll isional-radiati ve  coupling  in  the  negative 
glow  of  21 P  population  to  the  2^S  and  indirectly  to  the  2^S.  Some  of  the  21? 
production  results  in  metastable  production.  The  disagreement  in  the  highly 
abnormal  glow  discharge  of  1  .50  mA/cm‘L  is  in  part  due  to  the  assumption  of  no 
electron  reflection  from  the  anode.  The  Monte  Carlo  results  summarized  in 
Table  IV  show  a  very  high  average  kinetic  energy  (>100  eV)  for  electrons 
hitting  the  anode  in  the  highly  abnormal  glow  discharge.  Some  of  these 
electrons  should  be  reflected  at  the  anode.  The  reflection  from  the  anode  will 
increase  excitation  and  ionization  in  the  negative  glow  and  improve  agreement 
on  metastable  production  between  the  experiment  and  Monte  Carlo  simulations. 

The  reflection  will  not  significantly  change  the  current  balance  because  few  of 
the  reflected  electrons  will  return  to  the  cathode  fall  region  and  produce 
ionization . 

Figure  13c  shows  typically  3  metastable  atoms  produced  per  net  electron 
emitted  from  the  cathode.  This  number  can  be  used  to  estimate  the  importance 
of  metastable  bombardment  in  electron  emission  from  a  cold  cathode.  If  50$  of 


the  metastables  diffuse  back  to  the  cathode,  and  if  the  electron  emission 
coefficient  is  0.5,  then  the  metastables  could  be  responsible  for  75$  of  the 
electron  emission  from  the  cathode.^ 

The  overall  power  balance  in  the  discharge  is  of  considerable  interest. 
The  power  input  to  the  discharge  is  very  accurately  determined  from  the 
discharge  voltage  Vym  and  current  density  J^.  The  power  is  dissipated 
primarily  by  ion  transport  in  the  cathode  fall  region.  Ionization 
( 25  eV/event)  and  excitation  (-21  eV/event)  are  significant  but  less  important 
than  ion  transport.  The  other  major  power  dissipation  mechanism  is  due  to 
energetic  electrons  hitting  the  anode. 

The  amount  of  power  dissipated  by  each  of  the  four  mechanisms  is 
determined  using  the  Monte  Carlo  results  and  compared  to  the  input  power .  The 
Monte  Carlo  results  based  on  the  anisotropic  angular  distribution  for  ionizing 
collisions  are  used.  The  empirical  total  discharge  current  density  must  also 
be  used  in  analyzing  the  power  balance  because  the  Monte  Carlo  results  are  all 
per  (net)  electron  emitted  from  the  cathode.  A  simple  analytic  model  of  ion 
transport  in  the  cathode  fall  region  shows  that  it  consumes  very  nearly  60$  of 
the  power  dissipated.  Let  be  the  power  per  unit  area  dissipated  by  ion 
transport.  The  linear  variation  of  electric  field  and  the  equilibrium  charge 
exchange  drift  velocity  indicate  that 

d  d  _ 

Hit  -  I0'  J+  E  dz  =  /  c  J°/(d  -z)/d  E°  (d  -z )/d  dz  = 
o  o  ^ 

dc  o.u  =  J°  Vvm  0.8  =  JDVvm(0.8  J?/JD)  («U 

The  ratio  of  the  ion  current  density  at  the  cathode  J°  to  the  total  discharge 
current  density  JD  is  approximately  0.75,  this  ratio  is  simply  related  to  the 
current  balance  at  the  cathode  J°/(JD-J°)  as  determined  in  the  Monte  Carlo 
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simulations  .  Nearly  601  of  the  input  power  is  dissipated  by  ion  transport  at 
all  discharge  current  densities. 

The  power  per  unit  area  dissipated  in  ionization  (25  eV/event)  and 
excitation  (21  eV/event)  is  proportional  to  the  electron  current  density  at  the 
cathode  and  the  number  of  direct  ionization  and  total  excitation  events  per 
(net)  electron  emitted  from  the  cathode.  Associative  ionization  is  here 
included  with  total  singlet  plus  triplet  excitation.  The  expressions  are 

H ion  =  (25eV) (JD/e)[ ( JD~J°)/JD]  (direct  ionization  events/electron)  (42) 

and 

Hexc  =  (21 eV ) (JQ/e)[ (JD~J°)/JD]  (total  excitation  events/electron) .  (43) 

The  ratio  of  electron  current  density  at  the  cathode  (JD~J°)  to  the  discharge 
current  density  is  simply  related  to  the  current  balance  as  determined  in 
the  Monte  Carlo  simulations. 

The  power  dissipated  in  the  anode  per  unit  area  is 

^an  =  ( J[/e)[ ( total  ionization  events/electron)<T  >  (44) 

where  <Tan>  is  the  average  kinetic  energy  of  electrons  striking  the  anode. 

Again  all  of  the  quantities  are  from  the  Monte  Carlo  simulations  using  the 
anisotropic  angular  distribution  for  ionizing  collisions  except  for  the 
empirical  total  discharge  current  density  Jq.  The  results  of  the  power  balance 
comparison  are  summarized  in  Table  V.  The  power  balance  is  quantitatively 
understood.  The  second  to  last  row  is  the  total  power  dissipated  per  unit  area 
as  determined  from  the  Monte  Carlo  simulations.  The  last  row  is  the  empirical 


input  power  per  unit  area.  The  excellent  agreement  between  the  last  two  rows 
can  be  interpreted  as  additional  confirmation  that  the  current  balance  is 


correct . 

The  considerable  and  detailed  experimental  and  theoretical  information  on 
the  He  discharge  system  described  in  this  work  enables  us  to  analyze  the 
negative  glow.  The  negative  glow  region  starts  a  distance  dc  from  the  cathode 
where  the  cathode  fall  field  extrapolates  to  zero.  It  extends  to  the  anode  in 
our  discharge  cell.  The  separation  of  the  electrodes  was  chosen  so  that  no 
positive  column  would  form.  The  negative  glow  region  is  characterized  by 
extremely  small  electric  fields  and  a  high  density  ne  of  low  energy  electrons 
at  temperature  Tg.  These  low  energy  electrons  are  not  described  in  the  Monte 
Carlo  simulations.  Results  from  the  Monte  Carlo  simulations  such  as  <Tan>,  the 
average  energy  of  electrons  striking  the  anode,  are  for  the  high  energy  beam 
electrons.  The  low  energy  electrons  are,  as  we  shall  show,  trapped  in  a 
potential  energy  well.  Monte  Carlo  simulations  are  not  well  suited  for 
describing  trapped  particles. 

The  analysis  of  each  empirical  singlet  metastable  map  as  described  in 
Section  VI  produced  an  effective  rate  <ne  for  metastable  spin  conversion  due  to 
low  energy  electrons  in  the  negative  glow.  This  process  was  discovered  by 
Phelps  and  the  rate  constant  was  measured  for  room  temperature 
(kgTg  =  0.025  eV )  electrons. ^  The  dependence  of  this  rate  constant  on  electron 
temperature  is  (Te)  because  the  cross  section  scales  roughly  as  the  inverse 
of  the  electron  energy. ^  Although  we  do  not  a  priori  know  the  temperature  of 
the  low  energy  electrons  in  the  negative  glow,  it  is  relatively  easy  to  put 
some  limits  on  it.  Surely  the  temperature  of  the  low  energy  electrons  must  be 
greater  than  gas  temperature.  An  upper  limit  on  the  electron  temperature  la 
determined  by  noting  that  the  exothermic  metastable  spin  conversion  reaction 


depletes  the  singlets,  but  only  until  the  relative  singlet  versus  triplet 
metastable  density  comes  into  thermal  equilibrium  with  the  low  energy  electron 
gas.  Thus  the  empirical  ratio  Ms/Mt  in  the  negative  glow  must  be  greater  than 
or  equal  to  (1 /3)exp(-0.79/kBTe) .  These  limits  on  Te  and  knowledge  of  <ne 
determine  a  range  of  possible  electron  temperatures  and  a  corresponding 
electron  density  for  each  temperature  as  given  by  the  expression 

<ne  =  ne[<m/0.025eV/kBTe  -  (Mt/3MS) 

<m/0.025eV/kBTe  exp(-0.79eV/kBTe) ] .  (45  ) 

This  expression  is  represented  by  the  bold  solid  lines  of  Fig.  14. 

The  rather  high  densities  of  low  energy  electrons  suggests  the  existence 
of  a  potential  energy  well  for  electrons,  and  suggests  a  simple  model  for  the 
negative  glow.  (Without  a  potential  energy  well  the  diffusion  current  of 
electrons  to  the  anode  would  be  by  far  too  large.)  The  negative  glow  in  these 
experiments  can  be  thought  of  as  a  plasma  in  which  ion-electron  pairs  are 
produced  by  beam  electrons  from  the  cathode  fall  region  and  lost  by  a  process 
analogous  to  ambipolar  diffusion  to  the  anode.  The  analogy  to  ambipolar 
diffusion  is  not  perfect  because  the  electron  flux  density  on  the  anode  must  be 
larger  than  the  ion  flux  density  on  the  anode  by  JB/e.  The  loss  rate  of  low 
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energy  electrons  should  be  approximately  (1/2)  p+  kBTp(n/2dng)' /c  where  p,  is 

the  low  field  ion  mobility  and  dng  is  the  thickness  of  the  negative  glow.  The 

factor  of  1/2  is  due  to  the  fact  that  the  negative  glow  is  assymmetric,  the  low 

energy  electrons  cannot  diffuse  against  the  cathode  field.  The  rate  in  the 

negative  glow  of  ionization  per  unit  volume  <P>  can  be  determined  from  the 

Monte  Carlo  simulations  of  Table  IV.  The  rate  <P>  is  a  spatially  averaged 
d  +d 

quantity  (/dc  n8  P  dz)/dng. 


Inclusion  of  some  electron  reflection  from  the 


anode  will  increase  <P>,  but  at  most  by  a  factor  of  two  and  only  in  the  highly 
abnormal  glow  discharge  at  1.50  mA/cm'1'.  This  model  provides  a  second  relation 
between  the  density  ne  and  temperature  Tg  of  the  low  energy  electrons  in  the 
negative  glow 

<P>  -  nei  kB  Te  (TT/2dng)2/e  .  (46) 

This  expression  is  represented  by  the  light  solid  lines  of  Fig.  14.  The 
intersection  of  the  light  and  bold  lines  determines  the  density  and  temperature 
of  the  low  energy  electrons  in  the  negative  glow. 

This  simple  model  neglects  recombination.  We  do  observe  He2  molecular 
emission  bands  from  the  negative  glow.  These  bands  are  conclusive  evidence 
that  some  recombination  occurs.  The  inclusion  of  recombination  in  Eq.  (46) 
will  modify  the  low  temperature,  high  density  portion  of  the  balance  equation 
which  is  plotted  as  a  light  solid  line  in  Fig.  14.  Recombination  will  not 
significantly  change  the  intersection  point.  The  importance  of  recombination 
is  assessed  by  substituting  the  gas  density  and  the  density  and  temperature  of 
the  low  energy  electrons  as  determined  from  the  simple  model  into  the  known  He 
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recombination  rate.  °  Recombination  is  at  most  4$  of  the  diffusion  loss  rate 
(1/2)y  +  kBTe(TT/2dng)2/e. 

The  occurrance  of  the  intersection  on  the  steep  portion  of  the  plots  of 
Eq.  45  is  an  indication  that  the  metastable  spin  conversion  reaction  1 6a  is 
proceeding  to  equilibrium.  The  relative  density  of  the  singlet  and  triplet 
metastable  is  in  Boltzmann  equilibrium  with  the  low  energy,  high  density 
electrons  in  the  negative  glow.  Under  such  conditions  the  empirical  ratio  of 
singlet  to  triplet  metastables  is  providing  the  temperature  of  the  negative 
glow  electrons.  The  accuracy  of  this  determination  of  the  electron  temperature 
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in  the  negative  glow  is  not  easy  to  assess.  We  estimate  that  the  uncertainty 
of  the  temperature  and  density  of  the  negative  glow  electrons  determined  in 
Fig.  14  is  no  worse  than  ±50?.  This  completes  the  model  of  the  negative  glow. 

The  justification  for  ignoring  ionization  produced  in  the  negative  glow 
when  using  the  Monte  Carlo  simulations  to  determine  the  current  balance  at  the 
cathode  is  now  obvious  The  field  direction  reverses  a  distance  dc  from  the 
cathode  at  the  boundary  between  the  cathode  fall  and  negative  glow  regions. 

The  field  has  a  magnitude  on  the  order  of  one  V/cm  in  the  negative  glow.  Ions 
produced  in  the  negative  glow  drift  in  a  process  analogous  to  ambipolar 
diffusion  toward  the  anode  in  the  weak  field  of  the  negative  glow. 

The  remaining  issue  is  how  the  ion  density  of  a  few  x109cm~3  in  the 
cathode  fall  joins  smoothly  to  the  ion  density  of  >10ncnf3  in  the  negative 
glow.  One  at  first  might  suggest  that  there  will  be  a  large  diffusion  of  ions 
toward  the  cathode  across  the  cathode  fall-negative  glow  boundary.  This  does 
not  occur.  Furthermore  diffusion  is  not  a  valid  concept  within  a  few  ion 
mean-free- paths  of  the  cathode  fall-negative  glow  boundary.  Diffusion  is  an 
intrinsically  hydrodynamic  concept;  it  implies  a  density  gradient  expansion  of 
a  continuity  equation.  The  region,  within  six  mean-free-paths  for  symmetric 
charge  exchange,  from  the  cathode  fall-negative  glow  boundary  is  a 
non-equilibrium  region  even  for  ions  as  shown  in  Section  IV  of  this  paper.  The 
solutions  to  the  Boltzmann  equation  for  ions  given  in  Section  IV  have  a 
diverging  ion  density  at  the  boundary.  This  divergence  is  a  result  of  the  cold 
gas  approximation  for  the  ion  production  term  and  for  charge  exchange  collision 
term  used  in  the  Boltzmann  equation.  The  divergence  is  unphysical,  but  it  is 
essential  to  note  that  large  ion  density  gradients  at  the  cathode  fall-negative 
glow  boundary  are  built  in  to  the  exact  solution  to  the  Boltzmann  equation  for 
ions.  There  is  no  theoretical  difficulty  in  smoothly  joining  the  uniform  (few 


xlO^cm  3)  ion  density  of  the  cathode  fall  region  to  a  much  higher  (>1011cm  3) 


ion  density  of  the  negative  glow.  A  physically  reasonable  source  term  for 
electron  impact  ionization,  such  as  those  given  by  a  Monte  Carlo  simulation  for 
electrons  peaks  near  the  cathode  fall-negative  glow  boundary.  When  such  a 
source  term  is  combined  with  the  Green's  function  for  ions  as  in  Eq.  9  the 
sharp  ion  density  gradient  at  the  cathode  fall-negative  glow  boundary  appears 
naturally . 


IX.  Summary,  Conclusions,  and  Future  Work 
Various  laser  diagnostics  are  used  to  study  the  cathode  fall  and  negative 
glow  regions  of  a  He  glow  discharge.  The  water-cooled  A1  cathode  is  operated 
over  a  range  of  current  densities  from  a  near  normal  cathode  fall  of  173V  to  a 
highly  abnormal  cathode  fall  of  600V.  Optogalvanic  detection  of  Rydberg  atoms 
is  used  to  map  the  electric  field  in  the  cathode  fall  region.  The  Doppler 
width  of  a  transition  between  non-Rydberg  levels  is  used  to  measure  the  gas 
temperature.  The  gas  density  is  significantly  reduced  in  the  abnormal  glow 
discharge  due  to  ion  charge  exchange  heating.  Electric  field  maps  and  gas 
density  measurements  are  combined  with  an  analytic  treatment  of  ion  transport 
in  the  cathode  fall  region  to  determine  the  current  balance  or  ratio  of  ion  to 
electron  current  at  the  cathode.  This  current  balance  is  3-3  and  is  nearly 
independent  of  total  discharge  current.  The  densities  of  He  2^  and  2% 
metastables  are  mapped  using  laser-induced  fluorescence  and  laser  absorption 
spectroscopy  for  absolute  calibration.  An  analysis  of  the  metastable  maps 
provides  absolute  metastable  production  rates  in  the  cathode  fall  and  negative 
glow  regions  and  provides  a  constraint  on  the  density  and  temperature  of  the 


low  energy  electrons  in  the  negative  glow.  The  constraint  arises  from  our 
observation  of  a  suppression  of  the  2^  metastable  density  in  the  negative 


glow.  This  suppression  is  due  to  metastable  spin  conversion  by  low  energy 
electron  collisions  in  the  negative  glow. 

The  empirical  current  balance  at  the  cathode  surface  and  the  metastable 
production  rates  are  compared  to  theoretical  results  from  Monte  Carlo 
simulations  of  electron  avalanches.  A  Monte  Carlo  code,  based  on  the  null 
collision  technique  for  nonuniform  fields  as  developed  by  Boeuf  and  Marode,  is 
used  to  simulate  electron  avalanches.  The  agreement  between  the  empirical 
current  balance  and  the  result  of  Monte  Carlo  simulations  is  good.  The 
inclusion  of  an  anisotropic  angular  distribution  for  ionizing  collisions 
results  in  excellent  agreement,  even  in  the  highly  abnormal  cathode  fall  region 
where  E/N  reaches  3000  Td.  The  agreement  between  the  empirical  metastable 
production  rates  and  the  results  of  Monte  Carlo  simulations  is  good. 

The  high  density  of  low  energy  electror-,  in  the  negative  glow  suggests  a 
simple  model  for  the  negative  glow  and  suggests  the  existence  of  field  reversal 
at  the  boundary  between  the  cathode  fall  and  negative  glow  regions.  The 
empirical  constraint  on  the  density  and  temperature  of  the  low  energy  electrons 
in  the  negative  glow  is  combined  with  an  ionization  rate  predicted  by  Monte 
Carlo  simulations  to  determine  the  density  and  temperature  of  the  low  energy 
electrons . 

The  comparisons  of  the  empirical  results  and  results  from  Monte  Carlo 
simulations  lead  us  to  conclude  that  excitation  and  ionization,  as  embodied  in 
the  current  balance,  are  quantitatively  understood  in  the  cold  cathode  He 
discharge.  The  power  balance  is  also  quantitatively  understood.  We  conclude 
that  the  vast  majority  of  ions  striking  the  cathode  are  produced  in  the  cathode 
fall  region.  There  is  no  significant  "diffusion  like"  ion  current  at  the 
boundary  between  the  cathode  fall  and  negative  glow  regions.  Diffusion  is  not 
even  an  applicable  concept  because  of  the  lack  of  hydrodynamic  equilibrium 


within  a  few  ion  mean-free-paths  of  the  cathode  fall-negative  glow  boundary. 


Our  analytic  nonequilibrium  solutions  to  the  Boltzmann  equation  for  ions  in  the 
cathode  fall  region  correctly  describes  the  steep  increase  in  ion  density  at 
the  cathode  fall-negative  glow  boundary. 

The  conclusion  that  ions  reaching  the  cathode  surface  are  produced  in  the 
high  field  cathode  fall  region  is  generally  valid.  A  Bohm  presheath,  if  it 
existed,  could  preaccelerate  ions  before  they  reach  the  cathode  fall-negative 
glow  boundary.  Symmetric  charge  exchange  will  prevent  the  development  of  any 
structure  analogous  to  a  Bohm  presheath  in  the  negative  glow.  It  is  premature 
to  conclude  that  the  current  balance  of  3-3  is  typical  for  gases  other  than  He. 
We  can  conclude  that  the  experimental  and  theoretical  techniques  used 
in  this  investigation  are  broadly  applicable. 

A  quantitative  microscopic  picture  of  the  cathode  fall  and  negative  glow 
regions  is  emerging  but  major  challenges  remain.  A  major  experimental 
challenge  is  to  measure  the  relative  importance  of,  and  coefficients  for, 
electron  emission  from  the  cathode  due  to  ion,  metastable  atom,  and  VUV  photon 
bombardment.  The  measurements  should  be  done  in-situ  with  a  relatively 
nonperturbing  spectroscopic  diagnostic.  We  believe  that  the  spatial  and 
temporal  characteristics  of  carefully  chosen  optogalvanic  effects  will  provide 
this  information  on  electron  emission.  A  convincing  and  fully  self-consistent 
analysis  of  the  optogalvanic  effect  will  likely  be  difficult.  Another 
experimental  challenge  is  to  perform  independent  measurements  of  the  density 
and  temperature  of  the  low  energy  electrons  in  the  negative  glow  region. 

Again,  it  should  be  possible  to  make  the  measurements  with  a  relatively 
nonperturbing  spectroscopic  diagnostic.  On  the  theoretical  side  there  are  many 
challenges.  The  ultimate  goal  is  to  develop  relatively  simple,  efficient 
models  which  provide  fully  self-consistent  fields.  The  models  must  also  have 
good  microscopic  detail  and  absolute  accuracy. 
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Fig.  1  Schematic  of  the  experimental  apparatus  showing  all  three  detection 
methods:  optogal vanic ,  fluorescence,  and  absorption  detection. 

Fig.  2  (a)  Theoretical  Stark  map  for  the  n=11  singlet  m=0  levels  of  He. 

(b)  Theoretical  intensities  of  selected  Stark  components  as  a  function 
of  electric  field. 

Fig.  3  Optogalvanic  Stark  spectra  in  the  3-50  Torr  He  discharge  operating  at 
I^OmA/cm^.  The  distance  from  the  cathode  in  cm  at  which  each  spectra 
is  taken  is  indicated  at  the  right. 

Fig.  Electric  field  as  a  function  of  distance  from  the  cathode  for  five 

current  densities,  all  at  3.50  Torr.  The  lines  are  linear-least-square 
fits  to  the  data.  The  anode  corresponds  to  the  right-hand  side  of  the 
figure . 

Fig.  5  Experimentally  derived  gas  density  as  a  function  of  discharge  power 
measured  at  a  constant  pressure  of  3.50  Torr. 

Fig.  6  Integrals  of  the  Eq.  6  distribution  function  verses  the  number  of 

mean-free-paths  from  the  source  at  z  near  the  origin:  (a)  flux  density, 

o 

(b)  average  velocity,  (c)  ratio  of  average  velocity  to  local  equilibrium 
drift  velocity,  (d)  density. 

Fig.  7  Helium  2  S  metastable  density  as  a  function  of  distance  from  the  cathode 
for  five  discharge  current  densities  at  3*50  Torr. 

Fig.  8  Helium  ?^S  metastable  density  as  a  function  of  distance  from  the  cathode 
for  five  discharge  current  densities  at  3-50  Torr. 

Fig.  9  Helium  2  S  metastable  density  for  JD  -  1. 18mA /cm  .  The  points  are  the 
experimental  data  and  the  dashed  curve  is  the  calculated  fit  to  the 


Fig.  10  Helium  2^S  metastable  density  for  JD  -  1.l8mA/cm2.  The  points  are  the 
experimental  data  and  the  dashed  curve  is  the  calculated  fit  to  the 
data . 

Fig.  11  Monte  Carlo  histograms  showing  the  number  of  a)  excitation  events  and  b) 
ionization  events  per  (net)  electron  emitted  from  the  cathode  as  a 
function  of  distance  from  the  cathode.  The  histogram  for  total 
excitation  is  subdivided  into  21S,  2-^S  and  21  P  excitation.  The 
histogram  for  total  ionization  is  subdivided  into  associative  and  direct 
ionization . 

Fig.  12  Ratio  of  ion  current  to  electron  current  at  the  cathode  surface 

J?/(JD"J+)  a9  a  function  of  discharge  current  density.  The  points  are 
the  empirically  derived  values;  the  dashed  line  is  the  Monte  Carlo 
results  assuming  the  isotropic  scattering  for  ionizing  collisions;  the 
solid  line  is  the  Monte  Carlo  results  incorporating  the  anisotropic 
scattering  for  ionizing  collisions. 

Fig.  13  (a)  Helium  2^  production,  (b)  Helium  2^S  production  and  (c)  total 

metastable  production  as  a  function  of  discharge  current  density.  The 
production  values  presented  are  spatially  averaged  production  per  (net) 
electron  released  from  the  cathode.  The  points  are  empirical  results 
with  error  bars  reflecting  only  the  uncertainty  of  the  metastable 
density  measurements.  The  solid  lines  are  Monte  Carlo  results  using  the 
anisotropic  angular  distributions  for  ionizing  collisions. 

Fig.  14  Negative  glow  electron  density  versus  electron  temperature  for  five 

discharge  current  densities.  The  bold  line  is  the  relationship  from  the 
analysis  of  the  metastable  maps.  The  light  solid  line  is  from  the 
ambipolar-like  diffusion  of  the  electron  ion  pairs  in  the  negative  glow. 


The  intersection  of  the  two  lines  determines  the  density  and  temperature 


of  the  low  energy  electrons  in  the  negative  glow, 


i 


m 


Table  I.  Summary  of  the  analysis  of  the  electric  field  and  gas  density 
measurements  to  determine  the  current  balance  at  the  cathode  surface.  The 


dimensionless  ratio  J°/(JD~J°)  is  the  number  of  ions  produced  in  the  cathode  fall 
region  per  (net)  electron  emitted  from  the  cathode. 


(mA/cnd) 

0.190 

0.519 

0.846 

1 .18 

1  .50 

E°  (kV/cm) 

0.897 

1  .426 

1  .870 

2.395 

3.017 

dc  (cm) 

0.382 

0.301 

0.282 

0.300 

0.396 

Vvm 

0.173 

0.21 1 

0.261 

0.356 

0.600 

Vef  (kV) 

0.171 

0.215 

0.264 

0.359 

0.597 

N  ( 1 01 ^  cm  3 ) 

11.2 

10.8 

10.3 

9.48 

8.01 

p+  (10-10  C/cm3) 

2.08 

4 . 20 

5.88 

7.07 

6.74 

<m+v2/2>°  (eV) 

1  .66 

2.82 

3.89 

5.71 

8.92 

o  ( 1 0~16  cm2) 

- 

- 

23.3 

22.1 

21.1 

<vz>°  (105  cm/s) 

7.12 

9.28 

10.9 

13.2 

16.5 

(mA/cm2) 

CO 

•=T 

o 

0.390 

0.641 

0.933 

1.11 

J°/(JD-J?) 

3.52 

3.02 

3-13 

3.78 

2.85 
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Table  LI.  Summary  for  all  current  densities  of  parameters  used  in  the  model  as 
well  as  the  key  quantities  JPsdz,  Jptdz  and  <ne.  The  dimensionless  ratios 
ejPydz/' ( and  e/Ptdz/(JD-J°)  are  the  average  number  of  singlet  and  triplet 
metastables  produced  per  (net)  electron  emitted  from  the  cathode. 


j  (mA/cm2) 

0.190 

0.519 

0.846 

1  .18 

1  .50 

N  (10lD  cm~3) 

11.2 

10.8 

10.3 

9.48 

8.01 

Ds  (em2/s) 

127 

1  A0 

152 

173 

236 

Dt  (cmVs) 

136 

150 

163 

185 

253 

y  (10  1Cj  cm3/s) 

6.0 

6.5 

6.9 

7.8 

9.9 

B  (10  ^  cm'/s) 

3.5 

3.5 

3.5 

3.5 

3.5 

P3dz  ( 1  01  ^  crrf 2  s_1  ) 

0.454 

0.664 

0.830 

1 .05 

1 .38 

jptdz  (lO1*3  cm-2  s  1  ) 

0.71 1 

1  .86 

2.77 

3.56 

4.34 

eJ 

frsdz/(JD-J°) 

1.73 

0.824 

0.648 

0.680 

0,566 

GJ 

[ptdz/(JD-J°) 

2.71 

2.31 

2.16 

2.31 

1 .78 

<r 

ie  (10ii  s"1) 

1  .5 

8.5 

7.0 

5.5 

3.0 

T  ible  III.  Results  of  Monte  Carlo  simulations  of  electron  avalanches  based  on  the 


isotropic  angular  distribution  for  ionizing  collisions.  All  results  are  per  (net) 
electron  emitted  from  the  cathode.  The  ionization  in  the  cathode  fall  region 
includes  direct  plus  associative  ionization  occurring  between  the  surface  of  the 
cathode  and  a  distance  dc  from  the  cathode.  Other  quantities  include  all  events 
between  the  cathode  and  anode. 


(mA/cm‘ ) 

0.190 

0.519 

0.846 

1.18 

1  .50 

E°  (kV/cm) 

0.897 

1 .1426 

1 .870 

2.395 

3.017 

dQ  (cm) 

0.382 

0.301 

0.282 

0.300 

0.396 

N  (1016  cm'3) 

11  .2 

10.8 

10.3 

9.48 

8.01 

2^  production 

0.88 

0.88 

0.99 

1  .06 

0.96 

23S  production 

2.05 

1  .92 

2.04 

2.04 

1  .62 

21 P  production 

3.06 

3.37 

3.91 

4.50 

4.51 

Associative  ionization 

0.36 

0.38 

0.42 

0.47 

0.49 

Direct  ionization 

4.70 

5.98 

7.31 

9.02 

9.84 

Ionization  in  the 

cathode  fall  region 

3.89 

3.77 

4.06 

4.52 

5.71 

i 
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Table  IV.  Results 

of  Monte  Carlo  simulations  of 

electron  avalanches  based  on 

the 

anisotropic  angular 

distribution 

for  ionizing  collisions.  The 

quantity  <Tan> 

is 

the  average  kinetic 

energy  of  electrons  hitting 

the  anode.  Other  quantities  are 

per  (net)  electron 

emitted  from 

the  cathode. 

(mA/cm2) 

0.190 

0.519 

0.846 

1 .18 

1  .50 

E°  (kV/cm) 

0.897 

1  .426 

1 .870 

2.395 

3.0V 

dc  (cm) 

0.382 

0.301 

0.282 

0.300 

0. 39( 

N  (1016  cm-3) 

11.2 

10.8 

10.3 

9.48 

8.01 

21 S  production 

0.77 

0.75 

0.75 

0.74 

0.53 

23S  production 

1  .79 

1  .58 

1  .52 

1  .40 

0.83 

2^P  production 

2.82 

2.98 

3.14 

3.28 

2.51 

Associative  ionization  0.29 

0.31 

0.32 

0.34 

0.25 

Direct  ionization 

4.37 

5.37 

6.13 

6.88 

5.86 

<Tan>  (eV) 

14.0 

16.7 

21  .8 

36.5 

106.5 

Ionization  in  the 

cathode  fall  region 

3.40 

3.10 

2.97 

3.11 

3.24 

PI 


3 
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Table  V.  Power  balance  comparison  between  the  results  of  Monte  Carlo  simulations 


and  the  empirical 

H_  . 

■Tn 

input  power. 

The  variable  H  is 

the  sum  of  H^t, 

Hion*  Hexc’ 

and 

O 

Jp  (mA/cm  ) 

0.190 

0.519 

0.846 

1.18 

1  .50 

H^t  (mW/cm‘“) 

20.3 

66 

132 

254 

550 

Hion  (mW/cm2) 

4.7 

17 

33 

49 

52 

Hexc  (mW/cm2) 

5.1 

15 

26 

35 

31 

Han  (mW/cm2) 

3.4 

14 

35 

86 

268 

H  (mW/cm2) 

33.5 

112 

226 

424 

901 

VvmJD  (mW/cm2) 

32.9 

110 

221 

420 

900 

*»1.'  *  '■  ■  * 
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UfluojnjniMiLkiKt 


Schematic  of  the  experimental  apparatus  showing  all  three  detection 


(a)  Theoretics 


Calibration 

Fringes 


0.025 


0.076 


0.127 


0.178 


0.229 

0.280 

0.331 

0.382 


LASER  FREQUENCY  OFFSET  (cm"’) 


‘  ,T  r  ‘VV  ^  V 


.  f  .  ■  Vli  ■  ivivt* top jfv,y  v.y/-1  tswtp  v  r 

■  -  ■..»  * 1  *  ».**•*• ** yi tr y •. ,4  ■  ;*  « 


lorrespcnds  to  the  rignt-hand  side  of  the 


Density  (1016/cm3) 


Discharge  Power  (Watts) 


Fig.  •)  Experimentally  derived  gas  density  as  a  function  of  discharge  power 
measured  at  a  constant  pressure  of  3-50  Torr. 


Fig.  7  Helium  2’s  metastable  density  as  a  function  of  distance  from  the  cathode 
for  five  discharge  current  densities  at  3.50  Torn. 


discharge  current  densities  at  3.50  Torr. 


experimental  data  arid  fh°  dashpd  curve  13  the  calculated  fit  to  the 
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(a) 

^=0.519  mA/cm2f 


is®  associative  ionization 
□  direct  ionization 


■  ~  \\\W  »'  \WW'  W\\\'  \\\\\'  \W  \W\V  >W  WX'  >\W\\\\V  v\WV  WWW  'WWW,  .w\A? 


^=0.519  mA/cm2 


El  2’S  excitation 
E2  23S  excitation 
□  2’P  excitation 


Distance  From  Cathode  (cm) 


Fig.  11  Monte  Carlo  histograms  showing  the  number  of  a)  excitation  events  and  b) 
Ionization  events  per  (net)  electron  emitted  from  the  cathode  as  a 
function  of  distance  from  the  cathode.  The  histogram  for  total 
excitation  Is  subdivided  Info  7*3,  and  ?'p  excitation.  The 
histogram  for  total  Ionization  Is  subdivided  Into  associative  and  direct 
1 onl zat Ion . 


solid  line  Is  the  Monte  Carlo  results  Incorporating  the  anisotropic 


with  error  bars  reflecting  only  the  uncertainty  of  the  metastable 
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